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Abstract 

The authority of J. A. Wheeler in many areas of gravitational physics is immense, and there is a con- 
nection with the study of relic gravitational waves as well. I begin with a brief description of Wheeler's 
influence on this study. One part of the paper is essentially a detailed justification of the very existence of 
relic gravitational waves, account of their properties related to the quantum-mechanical origin, derivation of 
the expected magnitude of their effects, and reasoning why they should be detectable in the relatively near 
future. This line of argument includes the comparison of relic gravitational waves with density perturbations 
of quantum-mechanical origin, and the severe criticism of methods and predictions of inflationary theory. 
Another part of the paper is devoted to active searches for relic gravitational waves in cosmic microwave 
background radiation (CMB). Here, the emphasis is on the temperature-polarization TE cross-correlation 
function of CMB. The expected numerical level of the correlation, its sign, statistics, and the most appro- 
priate interval of angular scales are identified. Other correlation functions are also considered. The overall 
conclusion is such that the observational discovery of relic gravitational waves looks like the matter of a few 
coming years, rather than a few decades. 
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I. INTRODUCTION 



It is my honor and pleasure to be a lecturer at the first course of J. A. Wheeler School on 
Astrophysical Relativity. Wheeler is one of the founding fathers of the field of gravitational physics, 
and many of us are strongly influenced by his work and his personality. In particular, we often 
find inspiration in the great textbook of him and his colleagues 

I was fortunate to speak with J. A. Wheeler right before the publication of my first papers on 
relic gravitational waves Q]. That conversation helped me to shape my views on the subject. Being 
a young researcher, I was somewhat nervous about evaluation of my work by the towering scientist, 
but to my relief, Wheeler quickly understood the work and agreed with it. In what follows, we will 
be discussing relic gravitational waves systematically and in details, but I would like to start from 
describing my initial doubts and worries, and how Wheeler helped me to see them in a different 
light. 

I showed that the wave-equation for a gravitational wave h(7],x) = [/x(?7)/a(r/)]e*"''' in a homo- 
geneous isotropic universe with the scale factor a{r]) can be written as a Schrodinger-like equation 



It follows from this equation that the fate of the wave with the wavenumber n depends on the 
comparative values of n? and the effective potential U{7]) = a /a. If is much larger than |C/(?7)|, 
the wave does not feel the potential and propagates with the adiabatically changing amplitude 
h oc 1/a. In the opposite limit, the interaction with the potential is strong and the wave changes 
dramatically. The amplitude of the initial wave gets amplified over and above the adiabatic law 
h oc 1/a, and at the same time a wave propagating in the opposite direction is being created. 
This process results effectively in the production of standing waves. I called this phenomenon the 
superadiabatic amplification. Over the years, other cosmological wave equations were also modeled 
on Eq.([T]). In terms of the variable h{r]), where /i(r/) = fj,{rj) / a{rj) , Eq.([T|) has the form 



While analyzing various scale factors 0(77) and potentials U{ri), I was imagining them as being 
'drawn by a hand'. My first example was the potential U^rj) architypical for quantum mechanics 
- a rectangular barrier. This means that U{r]) = const in some interval of ry-time between rja and 
ijh, while U{ri) = outside this interval. I have shown that the waves interacting with this barrier 
are necessarily amplified. However, the postulated U {rj) caused some concerns. 




(1) 



h" + 2—h' + n^h = 0. 



(2) 



a 
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It is easy to make U{r]) = before rja and after rjb- Indeed, if a /a = 0, the scale factor 
0(77) is either a constant, hke in a flat Minkowski world, or is proportional to rj, like in a radiation- 
dominated universe. It is not difficult to imagine that the Universe was radiation-dominated before 
and after some crucial interval of evolution. However, if it is assumed that a" /a = const between 
rja and it is not easy to find a justification for this evolution, as the scale factor a(ri) should 
depend exponentially on ry-time in this interval. In terms of t-time, which is related to rj-time by 
c dt = a dr], the scale factor a{t) should be proportional to t. The Einstein equations allow this law 
of expansion, but they demand that the 'matter' driving this evolution should have the effective 
equation of state p = — (l/3)e. Some other potentials U (r/) drawn by a hand do also require strange 
equations of state. 

Having shown the inevitability of superadiabatic amplification, as soon as a" /a 7^ 0, I was 
somewhat embarassed by the fact that in some parts of my study I operated with scale factors 
drawn by a hand and driven by matter with unusual equations of state. Although equations of 
state with negative pressure had already been an element of cosmological research, notably in the 
work of A. D. Sakharov [3], I feared that Wheeler may dislike this idea and may say 'forget it'. 
To my surprise, he accepted the approach and even suggested a wonderful name: the 'engine- 
driven cosmology'. The implication was that although we may not know the nature of the 'engine' 
which drives a particular a{rj), this knowledge is not, for now, our high priority. Being inspired by 
Wheeler's attitude, I hurried to include the notion of the engine-driven cosmology, with reference 
to Wheeler, in the very first paper on the subject Q]. 

It is interesting to note that E. Schrodinger |J] felt uneasy about wave solutions in an expanding 
universe. (I became aware of his work much after the time of my first publications.) Schrodinger 
identified the crucial notion of the "mutual adulteration of positive and negative frequency terms 
in the course of time" . He was thinking about electromagnetic waves and he called the prospect 
of photon creation in an expanding universe an "alarming phenomenon" . From the position of our 
present knowledge, we can say that Schrodinger was right to be doubtful. Indeed, electromagnetic 
waves cannot be amplified and photons cannot be created in a nonstationary universe. Even though 
the wavelengths of electromagnetic and gravitational waves change in exactly the same manner, 
their interactions with external gravitational field are drastically different. The corresponding 
effective potential U{r]) in the Maxwell equations is strictly zero and the "alarming phenomenon" 
does not take place. All physical fields are tremendously 'stretched' by expansion, but only some 
of them are being amplified. 

In his paper, Schrodinger was operating with a variant of scalar electrodynamics, wherein the 
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coupling of scalar fields to gravity is ambiguous and can be chosen in such a way that the wave 
equation becomes identical to Eq. ([T]) , making the ainplification of scalar waves possible (for more 
details, see introductory part to Ref.P|). L. Parker {g] undertook a systematic study of the quan- 
tized version of test scalar fields in FLRW (Friedmann- Lemaitre-Robertson- Walker) cosmologies. 
As for the gravitational waves, there is no ambiguity in their coupling to gravity since the coupling 
follows directly from the Einstein equations. As we see, the "alarming phenomenon" does indeed 
take place for gravitational waves [2]. (The authors of publications preceeding Ref.^] explicitely 
denied the possibility of graviton creation in FLRW universes. In the end, it was only Ya. B. 
Zeldovich who wrote to me: "Thank you for your goal in my net".) 

It is important to realize that the possibility of generation of relic gravitational waves relies 
only on the validity of general relativity and quantum mechanics. The governing principles are 
part of the well-understood and tested physics. The underlying equation ([1]) admits an analogy 
with the Schrodinger equation (outlined above), but it can also be viewed as an equation for 
a classical oscillator with variable frequency. The phenomenon of superadiabatic (parametric) 
amplification of the waves' zero-point quantum oscillations is at the heart of the cosmological 
generating mechanism. In order to better appreciate this phenomenon, we shall briefly review the 
closely related laboratory-type problem of parametric amplification in a classical pendulum. 

Let us consider an ideal pendulum hanging in a constant gravitational field characterized by the 
free-fall acceleration g (see Fig(T]). The frequency of the oscillator is given by wq = ^/g/l, where I 
is the length of the pendulum, 

X + yX = 0. 

The amplitude of oscillations can be enhanced either by force acting directly on the pendulum's 
mass or by an influence which changes a parameter of the oscillator - in this case, its frequency. 
The equation for small horizontal displacements x{t) of the oscillator's mass takes the form 

x + uj'^{t)x = 0. (3) 

The simplest parametric intervention makes the length / time-dependent, as shown in FigH^. 
Then, uj'^{t) in Eq.® is w^(t) = {g — l)/l{t). Note that even if the gravitational accelaration 
g remains constant, it still gets modified by the acceleration term / arising due to the variation 
of length l{t). So, the correct io'^{t) differs from the naive expectation cj^(t) = g/l{t). (For more 
details about parametrically excited oscillators, see [^.) In general, both g{t) and l{t) are functions 
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FIG. 1: Parametric amplification, (a) Variation of the length of the pcndulmii. (b) Increased amplitude of 
oscillations. 



of time, in which case = [g{t) — l]/l{t) and Eq.dS]) reads 



X + X 



9{t) I 



= 0. 



(4) 



l{t) I 



If the external influence on the oscillator is very slow, i. e. ^ the ratio of the 

slowly changing energy E{t) and frequency uj{t), 



will remain constant. This ratio is called an adiabatic invariant The quantity N is also a 
'number of quanta' in a classical oscillator (see paper by Ya. B. Zeldovich, signed by a pseudonym, 
on how quantum mechanics helps understand classical mechanics [9]). 

On the other hand, if the oscillator was subject to some interval of appropriate parametric 
influence, the amplitude of oscillations and the number of quanta N will significantly increase, as 
shown in FiglTJ). To get a significant effect, the function uj{t) does not have to be periodic, but in 
its Fourier spectrum there should be enough power at frequencies around ujq. The final frequency 
does not need to differ from the initial ujq. 

One can notice the striking analogy between Eq.([T]) and Eq.Q. This analogy extends further 
if one goes over from the dispacement x{t) to the dimensionless angle variable (j){t) realated to x{t) 
by <j){t) = x{t)/l{t). (Compare with the cosmological relationship h{r]) = fi{ri) / a{r]) .) Then, Eq.dl]) 
takes the form 




(5) 



+ 2-(/.+ 




(j) = 



(6) 
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which is strikingly similar to Eq.(l2]). In cosmological equations, the analog of the ratio g{t)/l{t) is 



There are two lessons to be learned from this discussion. First, the necessary condition for a 
significant amplification of the wave is the availability of a regime where the characteristic time 
I a/a' I of variation of the external gravitational field (represented by the scale factor a{rj)) becomes 
comparable and much shorter than the wave period 27r/n, 

a'(r?) 



n <C 



(7) 



a(r/) 

Before and after this regime, the wave may be a high-frequency wave, that is, it may satisfy the 
condition n ^> \a' [t]) / a{r])\. But if during some interval of evolution the opposite condition ([Tj) is 
satisfied, and a" /a ^ 0, the wave will be siiperadiabatically amplified, regardless of whether the 
model universe is expanding or contracting 2|. 

Second, a classical pendulum should initially be in a state of oscillations - excited state - in 
order to have a chance to be amplified by an external 'pump' infiuence. Otherwise, if it is initially 
at rest, i.e. hanging stright down, the time-dependent change oil ox g will not excite the oscillator, 
and the energy of oscillations will remain zero. However, in the quantum world, even if the oscillator 
is in the state of lowest energy (ground, or vacuum, state) it will inevitably possess the zero-point 
quantum oscillations. One can think of these zero-point quantum oscillations as those that are 
being amplified by the external pump (assuming, of course, that the oscillator is coupled to the 
pump). This is where the quantum mechanics enters the picture. The initial ground state of 
the parametrically excited oscillator evolves into a multi-quantum state, and the mean number 
of quanta grows, when condition ([7]) is satisfied [ll|. There is no quantum state lower than the 
ground state, so the engine-driven cosmology will necessarily bring the properly coupled quantum 
oscillator into an excited state. 

The concept of the engine-driven cosmology remains perfectly adequate for the present-day 
research. We do not know what governed the scale factor of the very early Universe. It could be 
a lucky version of the scalar field (inflation) or something much more sophisticated and dictated 
by the 'theory of everything'. The lack of this knowledge is not important for the time being. 
By observing relic gravitational waves we may not be able to determine at once the nature of the 
cosmological 'engine', but we will be able to determine the behaviour of the early Universe's scale 
factor. Therefore, we will gain unique information about the 'birth' of the Universe and its very 
early dynamical evolution. 
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II. COSMOLOGICAL OSCILLATORS 



The physics of lab oratory- type oscihators have direct relevance to cosmological oscillators. In 
cosmology, we normally consider a universe filled with some matter and slightly perturbed in all 
constituents. It is convenient to write the perturbed metric of a flat FLRW universe in the form 

ds^ = -c^dt^ + a?{t){5ij + hij)dx'dx^ = a^{r]) [-dr]'^ + {6ij + hij)dx''dx^] . (8) 

The six functions hij {rj, x) can be expanded over spatial Fourier harmonics e^*"'"^, where n is a 

dimensionless time- independent wave- vector, 

+00 



Pij (n) hn iv) e'"''' c„ +Pij in)hn {v)e c„ 



(9) 



This representation requires some explanations. The mode functions hn (v) obey differential 
equations that follow from the perturbed Einstein equations (as an example, look at Eq.([T])). The 
wavelength of the mode n is given by A = 27ra/n, where the wavenumber n is n = IS 
convenient to take today's scale factor a{r]jj) to be equal to the 'size of the Universe', that is, a{r]ii) = 
21hi where Ih = c/Hq is today's Hubble radius and Hq = H{r]R) is today's Hubble parameter. 
Then, for a fixed moment of time, today, we can also write the laboratory-type expression for the 
same wavelenghth A: A = 2-11 /k, where k has the dimensionality of inverse length and is related to 
the dimensionless nhy k = n/2lH- The wave whose length is equal to today's 'size of the Universe' 
has n = 2tt. It is assumed that wavelengths can always be measured by unchangeable laboratory 
standards. 

By systematically writing n and distinguishing it from k we essentially follow the original mo- 
tivations of E. M. Lifshitz (see, for example. In contrast, some modern authors think that 
Lifshitz needs to be 'simplified' and 'modernized'. They use only one letter k and write k every- 
where where Lifshitz was writing n. This created quite a mess in wavenumbers and wavelengths. 
Among them you will see physical, non-physical, coordinate, comoving, proper, etc. 

Moving to complex Fourier coefficients Cn, we note that they are some particular numbers, if 
the left-hand-side (l.h.s) of Eq.Q is a deterministic, even if arbitrarily complicated, function. In 
the rigorous quantum-mechanical version of the theory, these coefficients will be promoted to the 
status of quantum- mechanical annihilation and creation operators Cn, acting on some quantum 
states. In the CMB applications, we will treat, for simplicity, c^, as random numbers taken 
from some probability distributions. The factor l/-v/2n in Eq.Q is a useful insertion inspired by 
quantum field theories. The normalization constant C will be discussed later. 



The gravitational field polarization tensors Pij (n) deserve special attention. As we shall see 
below, the polarization properties of the CMB radiation - our final destination ~ are intimately 
connected with the structure of these metric polarization tensors. Polarization of CMB and polar- 
ization of metric perturbations is not simply a coincidence in the usage of the word polarization. 

s 

The polarization tensors Pij (n) have different forms depending on whether the functions hij{rj,'K) 
represent gravitational waves, rotational perturbations, or density perturbations. Each class of 
these perturbations have two polarization states, so s = 1,2 for each of them. In what follows we 
will be considering gravitational waves and density perturbations. 

In the case of gravitational waves, two independent linear polarization states can be described 
by two real polarization tensors 

1 2 

Pij (n) = lilj - numj, Pij (n) = kmj + rmlj, (10) 

where spatial vectors (1, m, n/n) are unit and mutually orthogonal vectors. The polarization tensors 
(flOl) satisfy the conditions 

Pij6'^=0, kjn' = 0, iijP'^ = 26s's- (11) 

Two circular polarization states are described by 

L l/i 2\ _R l/i 2\ 

PiJ= ^ (Pij +i P^jj > Pv= ^ (Pij -i Pij) ■ (12) 

The left and right polarizations interchange under a coordinate reflection (altering the sign of or 
m*). In other words, gravitational waves can support a chirality, or 'handedness'. 
In the case of density perturbations, the polarization tensors are 

Pij=\lfi3^ P'^=-^^ + ^^'r (13) 

These polarization tensors satisfy the last of the conditions pi|) . The polarization tensors p3|) 
remain unchanged under coordinate mirror reflections, so density perturbations cannot support 
handedness. 

It is important to stress that from the position of general relativity, cosmological density pertur- 
bations represented by metric perturbations hij with polarization structure (jl3p . can be viewed as 
'scalar', or spin-0, gravitational waves. Indeed, although in the absense of matter, i.e. for T^i, = 0, 
the linearised Einstein equations admit spin-0 wave solutions, i.e. solutions with the structure (|13p . 
these solutions do not carry energy, do not affect the relative motion of test masses, and can be 
nullified by coordinate transformations. It is only spin-2 wave solutions, i.e. solutions with the 
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structure (jlOp . that carry energy, affect the relative motion of test masses, and cannot be removed 
by coordinate transformations. These spin-2 solutions are called gravitational waves. 

However, in cosmology, the spin-0 solutions survive and become non-trivial, as soon as metric 
perturbations with the structure (jl3p are supported by non-vanishing matter perturbations, that 
is, when STj^ ^ 0. These solutions, which unite gravitational field and matter perturbations, 
are called cosmological density perturbations. Therefore, cosmological density perturbations and 
cosmological gravitational waves, although separate from the point of view of algebraic classification 
of the tensor hij, are not entirely disconnected entities that should be treated by different theories. 
On the contrary, they should be viewed as parts of a common set of gravitational (metric) degrees 
of freedom. This can be regarded as a physical principle that will later guide our choice of initial 
conditions for density perturbations. 

One more comment is in order. In our presentation, we will be consistently using the class of 
synchronous coordinate systems ([8]), that is, we assume that /iqo = and /iQj = 0. We are not 
losing anything in terms of physics but we gain significantly in terms of technical simplicity and 
universality of our approach to gravitational waves and density perturbations. In principle, one 
can work in arbitrary coordinates, assuming that all components of metric perturbations, including 
/loo and /loi, are non-zero. This will not bring you any real advantages, but will complicate calcu- 
lations and can mislead you in issues of interpretation. Especially if you attempt to compare, say, 
gravitational waves descibed in synchronous coordinates with density perturbations described in 
the 'Newtonian-gauge' or other 'gauge' coordinates. Nevertheless, since various gauges and gauge- 
invariant formalisms are popular in contemporary literature, we will indicate, where necessary, how 
our formulas would be modified had we used arbitrary coordinates. 

III. QUANTIZATION OF GRAVITATIONAL WAVES 

Cosmological gravitational waves exist in the absence of matter perturbations, i.e. for (5T^ = 0. 

s s 

For each wavenumber n and polarization state s = 1,2 the mode functions hn {'n)^{'n) =l^n (v) 
satisfy the familiar equation ([T]). We assume that each of gravitational- wave oscillators n was 
initially, at some r] = ijq, in its ground state. Although certain oscillators could be somewhat 
excited without violating our perturbative assumptions, we do not see physical justification for 
such non-vacuum initial states. 

There is no such thing as the ground (vacuum) state without explicitely indicating the Hamilto- 
nian for which the state is the ground state. Shortly, we will explicitly write down the Hamiltonian 
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for gravitational waves (and, later, for density perturbations too). Specifying the Hamiltonian will 

also make more precise the concept of normalization of the initial mode functions to the zero-point 

quantum oscillations, or in other words the normalization to a 'half of the quantum in each mode'. 

It is important to remember, however, that, qualitatively, we already know the answer [2|]. The 

energy of a gravitational wave with wavelength Aq, contained in a volume (Ao/2)^, is equal to a 

half of the quantum, i.e. iV = 1/2 in Eq.(l5]), if the amplitude of the wave is at the level 

[gE n Ipi 
V c-^ ao Xq 



where Ipi = \J Gh/c^, uq = «(%) and Aq = 27r ag/n. Eq. (jl4p defines the initial vacuum spectrum 
of the gravitational wave (g.w.) amplitudes: hi{n) oc n. 

Shifting the initial time r]Q up to the boundary between the adiabatic and superadiabatic regimes 
at r] = r]i, we derive the estimate hi ~ lpi/\i ~ IpiHi/c. Then, we use the constancy of the metric 
amplitude h throughout the long- wavelength (superadiabatic) regime n <^ a' /a, Eq.Q, that is, the 
regime in which the wave is 'under the barrier a' /a\ The constancy of h follows from the constancy 
of the dominant (first) term in the approximate long- wavelength solution [2|] to the equations ([1]), 

^ = C, + cJ% (15) 
a J 

This allows us to write hj ~ /ij, where hj is the estimate of h at the end of the superadiabatic 
regime. After having emerged from 'under the barrier a'/a' the wave will again behave adiabatically. 

Using initial conditions ()14p and evolving classical mode functions through all the barriers and 
intervals of adiabatic evolution, one can derive today's metric amplitudes as a function of frequency, 
i.e. today's amplitude spectrum. For a typical engine-driven expanding cosmology shown in Figl^l 
and for its associated barrier a'/a shown in Fig|31 one expects to arrive at today's spectrum 
qualitatively shown in FigU] (for more details, see jl2|). 

The today's amplitudes in different parts of the spectrum are mostly determined by the thickness 
of the barrier a'/a, that is, by the duration of time that a given mode n spent under the barrier. 
All modes start with the initial Ni = 1/2, but if the mode n enters the barrier (i.e. satisfies 
the condition Aj ~ c/Hi) when the scale factor was ai{n), and leaves the barrier (i.e satisfies the 
condition Xj ~ c/Hf) when the scale factor was af{n), the final number of quanta in the mode 
will be Nf{n) = {aj/ai)^ = (Xf/Xi)'^ = {Hi/Hj)^. The very high-frequency modes with n above 
the tip of the barrier, n > ni, remain in the adibatic regime throughout their evolution, and their 
energy is being renormalised to zero. Therefore, the resulting amplitude spectrum quickly drops 
to zero at the high-frequency end. 
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FIG. 3: The barrier a'/a built from the the scale factor 0(77) of Figl2] versus crucial wavenumbers n defined 
by this barrier. 

The power-law index of today's spectrum h{n) in an arbitrary narrow interval of frequencies, 
say, between Ua and rib, as outlined in FigHl is expressible in terms of the initial vacuum spectrum 
hi{n) oc n and the power-law indeces /? of the scale factor 

a(r?) oc |r/|i+^ 

which approximates exact cosmological evolution in short intervals of time when the left and the 
right sides of the relevant portion of the barrier a'/a were formed. Today's spectral index of /i(n) 
in the discussed narrow interval of frequencies is given by 

h{n) oc n n^+l^' n'^^+f^f^ oc n^+^'-^z, (16) 

where /3j and (3f refer to the left and to the right slopes of the barrier, respectively. 
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FIG. 4: The characteristic gravitational-wave amplitude h(n) today as a function of frequency n. The size 
of the amplitude is mostly determined by the thickness of the barrier a' /a in the place, where it is traversed 
by the wave with a given wavenumber n. The barrier from Figl3]is shown at the top of the figure. 

The vacuum spectrum processed only on the left slope of the barrier, i.e. considered at times 
before processing on the right slope of the barrier, is called the primordial spectrum hp{n): 

hp{n)o^n^+^\ (17) 

For some historical and notational reasons, one and the same spectral index 2(2 + in the power 
spectrum h?'{n) of primordial metric perturbations is called nt in the case of gravitational waves 
and ri-s — 1 in the case of density perturbations (more details below, Eq. (j65p ). We will use the 
notation 2(2 + = n — 1. 

If the cosmological evolution a{r]) is such that the barrier a' /a is 'one-sided' and has only 
the right slope, the determination of the primordial spectrum requires additional considerations. 
The under-barrier amplification of the waves is still taking place ^, but the lack of the initial 
high-frequency regime makes the definition of the initial amplitudes somewhat ambiguous. 

To summarise, for a given cosmological evolution a{rj), one can qualitatively predict today's 
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piecewise amplitudes and slopes of the spectrum h{n) without doing any particularly detailed 
calculations. A barrier, whose shape is more complicated than the one shown in FigUl would have 
resulted in a more complicated shape of the generated spectrum. On the other hand, from the 
measured g.w. spectrum one can, in principle, reconstruct cosmological evolution a(r/) |l3l]. 



We will now turn to rigorous calculations based on quantum theory. The traditional, but not 
obligatory [l^, approach to quantum theory begins with a classical Lagrangian. The Lagrangian 
for a gravitational- wave oscillator of frequency n, and for each of two polarizations s = 1,2, has 
the form ( {isl and references there): 

2 ^ 



L 



gw 



1 / a 

2n \ao 



(h') 



212 

n h 



(18) 



where 



h 



327r3 



PI 



(19) 



This Lagrangian can be derived from the total Hilbert-Einstein quadratic action, where both gravity 
and matter parts of the action are taken into account. The classical equation of motion derivable 
from the Lagrangian (llSp in terms of the variable h is Eq.([2]), and the equation of motion derivable 
in terms of the variable /i, where h = /u/a, is Eq.([T]). 

The canonical pair of position q and momentum p for the oscillator can be taken as 



dL 



gw 



dh' 



i ( ^ 

n \ao 



h'. 



(20) 



Then, the classical Hamiltonian Hg^ = pq' — Lg^ reads 



gw 



n 
2 



-) f'+l- 

a / V ao 



(21) 



We now promote q and p to the status of quantum-mechanical operators and denote them 
by bold-face letters. Since we are interested in the initial conditions imposed in the early high- 
frequency regime, i. e. when n ^ a'/a, we can write the following asymptotic expressions for the 
operators: 



2T 



(22) 



' h a 
2 ao 



-ce 



-in{ri-rio) _|_ (,tg*n(f?-%) 



(23) 
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The commutation relationships for q, p operators, and for the annihilation and creation c, 
operators, are given by 



[q, p] = ih, 



c, c 



The asymptotic expression for the Hamiltonian Hgyj takes the form 



Hgw = h n c^c. 



(24) 



Obviously, the quantum state |0) satisfying the condition 



c|0) = 



is the state of the lowest energy, i.e. the ground (vacuum) state of the Hamiltonian (|24p . At 77 = ryo 
we get the relationships 



(Olq^lO) = (0|p2|0) = |, AqAp=|. 

The root-mean-square value of q in the vacuum state |0) is Qrms 
number with (I19D we derive 



hr 



An 



h/2. Combining this 



(25) 



which agrees with the qualitative estimate ([14]) . The adopted notations require C = ^/IGttIpi in the 
general expression ([9|). 

Now we will have to discuss exact Hamiltonians. The Hamiltonian built on the canonical pair 
/i, dLgw/dfi' manifestly illustrates the underlying pair creation process for gravitational waves. 
Specifically, we introduce 

Q = \ Q ^ = 1^ =Q Q 

V on riy/n Lpi oQ a 

and write the classical Hamiltonian in the form 



gw 



The associated annihilation and creation operators are 



The full quantum Hamiltonian can be written as 



H(77) = ncf c + (Tcf +a*c^, 



(26) 
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where coupling to the external field is given by the function cr(r/) = (i/2)(a'/a). 

The interaction with the external field can be neglected when n ^ a' /a, and then the first term 
in ()26p dominates and represents a free oscillator. The early-time asymptotic expressions for the 
Heisenberg operators, 



enter into formulas 

One can note that the notion of the barrier a' /a is convenient when one thinks of the problem 
in terms of the interaction Hamiltonian (j26p and the first-order Heisenberg equations of motion, 
whereas the barrier a" /a is more convenient when one thinks of the problem in terms of the 
second-order Schrodinger-like equation ([1]). 

The emerging correlation between the traveling modes n and — n, which leads to the production 
of standing waves referred to in the Introduction, is described by the 2-mode Hamiltonian (see 11 1 
and references there): 

H(?y) = ncJjCn + ncLnC_n + 2o-cJjcLn + 2o-*CnC_n- (27) 



This Hamiltonian can be viewed as the sum of two Hamiltonians (|26p . The Hamiltonian (j26p is 
called a 1-mode Hamiltonian. 

The defined quantum-mechanical operators and Hamiltonians, plus the assumption about a 
particular initial state of the field, fully determine dynamical evolution of the field, its statistical 
properties, correlation functions, and eventually the observational predictions for later times. 

IV. SQUEEZING AND POWER SPECTRUM 

The quantum- mechanical Schrodinger evolution transforms the initial vacuum state |On)|0_n) 
into a 2-mode squeezed vacuum state, which is equivalent to a pair of 1-mode squeezed vacuum 



states (see ll[ and references there). In the Heisenberg picture, the initial state of the field does 
not evolve, but the operators do, and their evolution is ultimately described by the mode functions 

s 

hn (r/)- 

It is the variance of the oscillator's phase that is being strongly diminished (squeezed), whereas 
the variance of the amplitude is being strongly increased. The Gaussian nature of the initial 
vacuum state is maintained in the course of the Schrodinger evolution, but the variances of phase 
and amplitude in the resulting squeezed vacuum quantum state are dramatically different. The 
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parameter of squeezing and the mean number of 



regime, and stop growing only at its end 



quanta are growing all the way up in the amplifying 



(The multi-quantum nature of the developing 
squeezed vacuum quantum state is behind the continuing debate over the 'quantum-to-classical' 
transition, 'decoherence', etc.) The phenomenon of squeezing allows us to treat the resulting 
quantum states as a stochastic collection of standing waves. The squeezing and the associated 
picture of standing waves is very i.p„.an. „B«™ti„nan. fl. I. lead, .o o»I,a.„,. fea.u.. 
in the metric power spectrum and, as a consequence, to oscillatory features in the angular power 
spectrum of CMB temperature and polarization anisotropics. We will discuss these oscillations 
later. 

Having accepted the initial vacuum state |0) of the gravitational (metric) field ([9]), we can 
calculate the variance of the field: 



(0|/i,,(r/,x)/i*^(r/,x)|0) 



oo 

t^Jn^^lL M^'i- (28) 



The quantity 

( nn iri\ — 

27r2 



h\n,^) = ^n' J] |L(r/)P (29) 

s=l,2 



gives the mean-square value of the gravitational field perturbations in a logarithmic interval of n 
and is called the metric power spectrum. The spectrum of the root-mean-square (rms) amplitude 
h{n,r]) is determined by the square root of Eq. (i29]) . 

s 

Having evolved the classical mode functions hn {f]) some arbitrary instant of time rj one 

can find the spectrum h{n, rf) at that instant of time. As mentioned before, the spectrum calculated 
at times when the waves of today's interest were in their long-wavelength regime (that is, longer 
than the Hubble radius at that times) is called the primordial spectrum. The today's spectrum, 
i.e. spectrum calculated at = r//j, is normally expressed in terms of frequency v measured in Hz, 
v = uHq/Att = nvn/^T^- The spectral rms-amplitude is denoted by h^msi^)-, or simply h{i'). The 
mean-square value of the gravitational-wave field in some interval of frequencies between vi and 
U2 is given by the integral: 

1/2 



J h\u)^. (30) 



The spectral function h'^{u) depends on frequency ly, but is dimensionless. The dimensionality of 
Hz""*^ is carried by the function h?[v)/v. 
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For gravitational waves which are comfortably shorter than the Hubble radius, one can also 
calculate the spectral gravitational- wave energy density /Og^(z^)c^ and the total g.w. energy density 
in some interval of frequencies. These quantities are expressible in terms of h?{y): 

Pgw{y) = ^h?'{y)v'^, Pgw{j^l,l'2) = J Pgw{T^) — - (31) 

For the purpose of comparing a g.w. background with other sorts of matter, it is convenient to 
introduce the cosmological r2gt„-parameter and its spectral value ngw{i'). As with all other sorts 
of radiation, this parameter is defined by 

= (32) 

Pcrit 



where pcrit = '^Hq/^iItG. Using Eq. ([3T]) one can also write [1 



f!3,.(>') = I — I ^ (33) 



The total between frequencies vi, is given by 

^gw{v\,V2)= \ ^gw{v) . 

One should be weary of the confusing definition 

^gwiy) 



1 dpgy,{y) 



Pcrit 

often floating in the literature. As it stands, this relationship is incorrect. It can be made consistent 
with the correct definition ()32p only if one assumes that what is being differentiated in this formula 
is not the spectral density Pgwi^^), but a logarithmic integral of this quantity in the limits between 
some fixed i^i and a running v. 

To make more precise the expected qualitative graph for hrmsin-) in FigjH as well as previous 
theoretical graphs in Fig. 4 of Ref . {igI . we need to use some available observational data. We make 
the fundamental assumption that the observed CMB anisotropics are indeed caused by cosmological 
perturbations of quantum-mechanical origin. If so, the contribution of relic gravitational waves to 
the large-scale anisotropics should be of the same order of magnitude as the contribution of density 
perturbations (we will show this in more detail below). This allows us to determine the position 
and the slope of the function hrms{i^) at frequencies near the Hubble frequency z^// ~ 2 x 10^^^ Hz 
and then extrapolate the spectrum to higher frequencies. 

We choose hrmsii^n) and primordial spectral index n in such a way that the lower-order CMB 
multipoles produced by relic gravitational waves are at the level of the actually observed values. 
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Then, today's spectra for hrms{i^) with spectral indeces n = 1 (/3j = —2) and n = 1.2 = 
— 1.9) are shown in Figl2] (for more details, see Q and Q). At frequencies around i'h we have 



hrmsii^H) ~ 10^^, so that the present-day mean number of quanta Nf exceeds 10^°*^. 

When extrapolating the spectrum to higher frequencies we assume that the entire left slope of 
the barrier in Figl3]was formed by a single power-law scale factor with one and the same /3j = (3: 

a{r,) = loH^^^^ . (34) 

This seems to be a reasonable assumption given the relative featurelessness of the initial stage - 
the overall energy density was 10 orders of magnitude lower than the Planckian density and was 
barely changing. Specifically, at the i-stage, we consider two examples: (3 = —2 and (3 = —1.9. 
The power-law indeces 13 f at the matter-dominated and radiation-dominated stages are well known: 
f3f = 1 and /3j = 0, respectively. The effective pressure p and the energy density e of matter driving 
the general power-law evolution ()34p with a given constant /3 are related by the effective equation 
of state p = [(1 - /3)/3(l + f3)]e. 

As for the high-frequency part of the spectrum, it was calculated under the assumption that 
the z-stage of evolution shown in Figl2] was governed by matter with a stiff equation of state 
p = . = -(1/2)) advocated long ago by Zeldovi*. The i..e of the bao. .action of the 
created gravitons on the "pump" field a{7]) becomes important for this sort of values of /3j [0], 



18l |. Of course, the existence of such an interval in the past evolution of the very early Universe 
cannot be guaranteed. In any case, the waves with frequencies above 10^*^ Hz have never been in 
superadiabatic regime, they remain in the vacuum state. The renormalization (subtraction of a 
"half of the quantum" from each mode) cuts off the spectrum at these high frequencies. At lower 
frequencies, the renormalization has practically no effect on the spectrum. 

The spectra of ilgw{i'), shown in FiglHl are derived from hrms{t^) according to Eg. ([33]) . As 
was already mentioned, the substantial rise of the spectrum at very high frequencies cannot be 
guaranteed, and the shown graphs in this area of frequencies should be regarded as the upper 
allowed limits (especially for the model with the primordial spectral index n = 1.2). 

There are two comments to be made about Figl5]and Figl6j First, the phenomenon of squeezing 
and standing waves production is reflected in the oscillations of the metric power spectrum and 
0-spectrum as functions of frequency v. The first few cycles of these oscillations are shown in 
the graph of hrms{i^)- In the CMB sections of the paper, we will show how these oscillations that 
existed in the recombination era translate into peaks and dips of the angular power spectra for 
CMB temperature and polarization anisotropies observed today. 
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Frequencv V, Hz 



FIG. 5: The present-day spectrum for hrms{i^)- The sohd hne corresponds to the primordial spectral index 
/3 = —1.9, i.e. n = 1.2, while the dashed line is for /3 = —2, i.e. n = 1. 

Second, since the primordial spectrum of quantum-mechanically generated perturbations has 
the form (see Eq. (fT7|) l: 

/i2(n)ocn2(/^+2)^ (35) 

theoretical considerations suggest that the preferred range of values for the primordial spectral 
index n is n > 1, i.e. preferred spectra are 'blue'. Indeed, for the opposite range n < 1, the 
primordial spectra would be 'red', and the integral (j28p would be power-law divergent at the 
lower limit. To avoid the infra-red divergency one would need to make extra assumptions about 
bending the spectrum down in the region of very small wavenumbers n. In its turn, a primordial 
spectral index n > 1 requires that the power-law index (3 of the initial scale factor a(r/) should be 
/? > —2. This is because j3 and n are related by n = 2/3 -|- 5. The value n = 1 (/3 = —2) describes 
the spectrum known as a flat, or Harrison-Zeldovich-Peebles, or scale-invariant spectrum. It was 
originally proposed in the context of density perturbations. The spectral index n = 1 marks the 
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Frequency v, Hz 



FIG. 6: The present-day spectrum for flgwlv). The solid hne corresponds to the primordial spectral index 
/3 = —1.9, i.e. n = 1.2, while the dashed line is for /3 = —2, i.e. n = 1. 

beginning of the trouble ~ the metric variance (|28p is logarithmically divergent at the lower limit. 

V. DENSITY PERTURBATIONS 

To discover relic gravitational waves in CMB anisotropics we have to distinguish their effects 
from other possible sources of CMB anisotropics, and first of all from density perturbations. In 
comparison with gravitational waves, the possibility of quantum-mechanical generation of density 
perturbations requires an extra hypothesis, namely, the appropriate parametric coupling of the 
matter field to external gravity. But we can assume that this hypothesis was satisfied by, for 
example, a lucky version of the scalar field. Can the amplitudes of quantum-mechanically generated 
density perturbations be many orders of magnitude larger than the amplitudes of relic gravitational 
waves ? This would make the search for relic gravitational waves practically hopeless. We will 
show that this can never happen. Density perturbations can compete with relic gravitational waves 



20 



in producing the large-scale CMB anisotropies, but can never be overwhelming. 

We are mostly interested in the gravitational field (metric) sector of density perturbations. 
Indeed, it is perturbations in the gravitational field that survived numerous transformations of the 
matter content of the Universe. Primordial matter (for example, scalar field) has decayed, together 
with its own perturbations, long ago. But it is primordial metric perturbations that have been 
inhereted by density perturbations at radiation-dominated and matter-dominated stages. 

As a driving 'engine' for the very early Universe, it is common to consider the so-called 
minimally-coupled scalar field (p{ri,x), with the energy-momentum tensor 



(36) 



The polarization structure p3|) of density perturbations allows us to write the n-mode metric 
perturbation: 

hij = h{7j)Q6ij + hi{r])n~'^Q^ij, (37) 

where Q = e^*"'^. Two polarization amplitudes h{r]) and /i;(ry) are accompanied by the third 
unknown function - the amplitude 9?i(??) of the scalar field perturbation: 

(p{r],x) = <fo{ri) +'Pi{v)Q- 

The perturbed Einstein equations allow us to find all three unknown functions. The important 
fact is that, for any V{ip), there exists only one second-order differential equation, of the same 
structure as Eq.([T]), that needs to be solved jl^: 



//^ 12 (0^/7)' 



0. (38) 



This equation coincides with the gravitational wave equation ([T]) if one makes there the replacement 



i{r]) a(7/)v^7(??). (39) 



We will call function /u satisfying Eg. ([38]) /ig, and that satisfying Eq.([T|) fix- 
The function 7(7/) is defined by 

For power- law scale factors (i34|) . this function reduces to a set of constants, 

^ 2 + /3 
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The constant 7 degenerates to zero in the hmit of the expansion law with /3 = —2, that is, in the 
Hmit of the gravitational pump field (an interval of deSitter evolution) which is responsible for the 
generation of cosmological perturbations with flat primordial spectrum n = 1. 
In terms of t-iime the function 7 is 



lit) 



H 



It is also related to the cosmological decelaration parameter q{t) = —aa/a'^: 7(t) = 1 + q{t). The 
unperturbed Einstein equations allow us to write 

,/\ 2 



(40) 



or, equivalently, 



H 



(41) 



where k = SvrG/c^. The function 7(t) is sometimes denoted in the literature by e(t). 

As soon at the appropriate solution to Eq. ()38p is found, all unknown functions are easily calcu- 
lable: 



h{rj) = -H{r,) 
c 



rv 



(42) 



a' 



h" - ^h' 
H' 



(43) 



V2k 



«V7 



h 



(44) 



The arbitrary constant Cj in Eq. (l42p reflects the remaining coordinate freedom within the class 
of synchronous coordinate systems ([8]). Indeed, a small coordinate transformation 



(45) 



generates a 'gauge transformation' 



hiT]) = hir]) + C^ = h{ri) + C 



(46) 



The term with constant C is automatically present in the properly written general solution ([l2 
for /i(ry). 
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By combining the transformation (j46p with its time-derivative, one can build a quantity which 
does not contain the constant C and therefore is a 'gauge-invariant' quantity with respect to 
transformations preserving synchronous coordinates. Specifically, one can check that 

h-h- = h-h-. 



Denoting this quantity by C, and taking into account Eq. (j42p . we write 

(,^h-h'^ = J^. (47) 

In terms of the variable C,{rj) the fundamental equation (|38p takes the form 

C" + 2^^^C' + n\ = 0. (48) 

Not surprisingly, this equation coincides with the gravitational- wave equation ([2]), if one replaces 
the gravitational-wave function h with C,, and a with a^. For density perturbations, the metric 
variable C, is the physically relevant quantity that plays the same role as the metric variable h for 
gravitational- wave perturbations. 

The fact that we are working with scalar gravitational waves supported by scalar field fluctu- 
ations, instead of normal tensor gravitational waves, has boiled down to the necessity of a single 
modification: the substitution (|39p in the gravitational-wave equations. 

It is appropriate to say a few words about cosmological gauge transformations in general. Their 
ori gin is related to the notion of Lie transport performed on a manifold covered by some coordinates 



[20|. Lie transport is being carried out along a given vector field .^'^(x"). This is a quite formal 
construction which respects only the transformation properties of fields defined on the manifold. 
These fields are not required to satisfy any physical equations. An infinitesimal Lie transport 
changes the field by an amount equal to its Lie derivative. This change can be viewed as a rule by 
which new values of the field are assigned to the same point x". It is only with some reservation 
that we borrow the name gauge transformation from physical field theories and apply it to this, 
always valid, mathematical procedure. 

In gravitational applications, the vector field E,^{x") is usually associated with an infinitesimal 
coordinate transformation = x^ —^^{x") [lo| . Then, for example, an infinitesimal increment of 
the metric tensor g^uix") is given by the Lie derivative of g^u- It can be written as 

A particular case of these transformations is represented by our Eqs. (j45^ I46 p . 

23 



Since we operate with only four components of an arbitrary vector and, potentially, with 
many fields transforming with the same one can build some combinations of these transformed 
fields, in which the functions cancel out. Therefore, these combinations become gauge- invariant 
quantities. An example is given by our Eq. (j47p . 

Certainly, it is an exaggaration to claim that "only gauge-invariant quantities have any inherent 
physical meaning". It is as if you were told that where specifically you are going in your car 
has no physical meaning, because the components of velocity are coordinate-dependent, and it is 
only the readings of your speedometer that have an inherent physical meaning, because they are 
coordinate-independent. In any case, there exists an infinite number of gauge-invariant quantities. 
For example, a product of a gauge-invariant quantity with any 'background' function produces 
a new gauge-invariant quantity. The fact that a quantity is gauge-invariant does not answer the 
question of its physical interpretation. Nevertheless, it is useful to know gauge-invariant quantities. 

Had we started from arbitrary coordinates, the perturbed metric dS]) would have contained two 
more unknown functions, namely A{r]) and B{r]), [2l| 

9r,r, = + 2A{i])Q), g^i = a\r])B{r])-Q^i. 

n 

The ry-transformation generalizing our Eq. (|45p would read 

fj = r] + T{r])Q, (49) 

where T(r/) is an arbitrary function of rj. The perturbed metric components h and A would 
transform as follows (any possible x*-transformations do not participate in these relationships): 

h = h- 2-r, A = A-T' --T. (50) 
a a 

The functions T and T' cancel out in the gauge-invariant metric combination 

Cg = h- —h . 



Obviously, Q reduces to C for transformations (I45p preserving synchronous conditions, that is, 
when A = A = 0. The previously introduced [22] quantity Cbst, where BST stands for Bardeen, 
Steinhardt, Turner, can also be reduced, after some work, to our up to a coefficient —(1/2). 
Since we assume that in addition to the tensor metric field g^u, a scalar field 

(/?(??, x) = ipo{ri) + 6ip{ri)Q 
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is also defined on the manifold, we can write down its gauge transformation under the action of 
Eq.dMD: 

I^ = 5^ + ^'oT. (51) 

The arbitrary function T(?y) cancels out in the gauge-invariant combination 

V{r,)=5^ + ^^^'oh. (52) 

Although the combination (j52p is often quoted in the literature, this object is something like a 
'half of a horse plus half of a cow'. It combines physically separate quantities - one from metric 
another from matter - and its physical interpretation is obscure. However, using the unperturbed 
Einstein equation (j40p and solution (j44p for 5ip = (pi one can show that V reduces to the gauge- 
invariant metric perturbation ( times a background-dependent factor: 



The quantity \'2KaV is called Chibisov, Lukash, Mukhanov, Sasaki 



23l425l| variable uclms- 



ucLMS = aVlC- (53) 

In our approach, this variable is simply the function fj,s satisfying Eq. (j38p . modeled on Eq.([T]). 
The factor -^/7 in Eq. (j53p will be a matter of great attention when we come to the discussion of 
quantization procedures. 

Returning to the function C(^)) one can notice that this metric amplitude is practically constant 
in the regime when is much smaller than the effective potential (a^y^)" /a,^. This behaviour 
is similar to the constancy of the gravitational wave amplitude h = /ir/a throughout the long- 
wavelength regime. Indeed, in full analogy with the long- wavelength solution ()15p . the general 
solution to Eq. (|38l) in this regime reads 



For usually considered expanding cosmologies, the term with constant C2 is decreasing, and there- 
fore the dominant solution is /i ~ Ci for gravitational waves and (" ~ Ci for density perturbations. 
The constancy of C allows one to easily estimate the value of metric perturbations at much later 
times, at the radiation-dominated and matter-dominated stages, as soon as one knows the initial 
value of Q. 

The constancy of C, (when one can neglect the term with C2) is sometimes called a conservation 
law. This association is incorrect. Genuine conservation laws reflect symmetries of the system, and 
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conserved quantities are constants independently of the initial conditions. For example, the energy 
of a free oscillator is a constant independently of initial positions and velocities. In our problem, 
the genuine conservation law for ( would look like an empty statement = 0. 



VI. QUANTIZATION OF DENSITY PERTURBATIONS 

A quantum system is defined by its Hamiltonian. Whether or not the Hamiltonian follows from 
some assumed classical Lagrangian should not be a question of major concern [l^. From this point 
of view, the quantization of density perturbations, similarly to the quantization of gravitational 
waves (j26p . is defined by the Hamiltonian [2] 

H(r/) = nJc + fjc"''^ + (T*c2. (54) 

The pair creation of scalar gravitational waves C by the external pump field is regulated by the 
coupling function (T(r/) = {i/2)[{a.y^y /(ay/^)]. However, a more traditional approach begins with 
a Lagrangian, and here too, the Hamiltonian (I54p can be derived from quadratic perturbation 
terms in the total Hilbert-Einstein action, where both gravity and matter (scalar field) are taken 
into account. 

The Lagrangian for an n-mode of density perturbations, after some transformations of the total 
Hilbert-Einstein quadratic action, can be written in the form [l^ 



1 / aV7 ' ' 



^'^^ 2n V«o^/7o 



{Cf-n'c'\, (55) 



where 

and ao, 70 are values of the functions a(r/), 7(7/) at rj = r]Q where the initial conditions are being 
set. The Euler-Lagrange equations derivable from this Lagrangian in terms of ^ are Eq. (j48p . and 
in terms of fis ~ Eq . (p8]) . 

Obviously, the Lagrangian ()55p for coincides with the gravitational wave Lagrangian (llSp for h 
after the replacement of the factor a/ao with the factor 0^/7/00-^/70. Starting from the Lagrangian 
(|55p and building on the canonical pair dLdp/dfx'g one can derive the Hamiltonian (j54p by 
doing exactly the same steps that have led us from the Lagrangian (|18p to the Hamiltonian (|26p . 

It should be noted that the total classical Lagrangian (j55p admits, as always, some freedom of 
modifications without affecting the Euler-Lagrange equations. In particular, the derived equations 
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of motion ()48p . (j38p will remain exactly the same, if one changes the Lagrangian (j55p to a new one 
by multiplying (|55p with a constant, for example with a constant 70: 



L 



dpinew) 



1 / 



70 



(C')'-n^C 



1 / a 



2n\aQ ^ 



(57) 



2ra V^oVTo, 

This new Lagrangian degenerates to zero in the limit 7 — )• 0. We will discuss later the subtleties 
in quantum theory that arise after such a modification of the Lagrangian. 

Similarly to gravitational waves, we impose initial conditions in the early regime of a free 
oscillator, that is, when n ^ (o-y/j)' /(o-y/^)- We choose the canonical pair 



P 



1 / a\/7 



and write the asymptotic expressions for the operators: 

2 L 



(58) 



(59) 



2 aov/70 



-ce 



-in(ri-rio) _|_ ^in(T]-r!o) 



with 



[q, p] = ih, 



c, c' 



Obviously, a quantum state satisfying the condition 



c|0) = 



(60) 



(61) 



(62) 



is the ground (vacuum) state of the Hamiltonian (I54p . Calculating the mean square values of q 
and its canonically conjugate momentum p, we find 



n 



(0|q^|0) = (0|p^|0) =-, AqAp = - 



Returning to C from C according to Eq. (j56j) . we find 

Crms = ((0|C'|0)) 



An 



(63) 



that is, exactly the same value as the initial amplitudes (j25p for each of two polarization components 
of gravitational waves. 

Extrapolating the initial time r/o up to the boundary between the adiabatic and superadiabatic 
regimes at r] = rji, we derive the evaluation Crms ~ Ipi/^i- This evaluation, plus the constancy of 
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the quantity ^ throughout the long-wavelength regime, is the foundation of the result according 
to which the final (at the end of the long-wavelength regime) metric amplitudes of gravitational 
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waves and density perturbations should be roughly equal to each other 

The primordial (^-spectrum has the same form as the primordial gravitational-wave spectrum 
in Eq.d 



In this approximation, the spectral indeces are equal: 

ns-l = nt = 2{/3 + 2) = n-l. (65) 



Scalar fields (|36p can support cosmological scale factors with /3 in the interval — 1 < /3, /3 < —2. 
The ratio of (j35p to (j64p is approximately 1 for all spectral indeces near and including /3 = —2 
(n=l). 

Having strictly defined the dynamical equations and the initial (quantum ground state) con- 
ditions, one can calculate from formula ()29p the exact power spectrum of metric perturbations 

s 

associated with density perturbations. In this formula, two polarization components hn iv) 



now determined by the conventions p^ . ([37|) . and the constant C is C = \/247r/p;. 

The employed approximations cannot guarantee that in the real Universe the coefficients in 
Eq. psp and Eq. ()64p should be exactly equal to each other. But there is no reason for them to be 
different by more than a numerical factor of order 1. Therefore, the lower order CMB multipoles, 
induced primarily by metric perturbations which are still in the long- wavelength regime, should 
be approximately at the equal numerical levels for, both, density perturbations and gravitational 
waves. 

As we shall now see the infiationary theory differs from the described quantum theory of density 
perturbations by many orders of magnitude, let alone numerical coefficients of order 1. The dis- 
crepancy becomes infinitely large in the limit of the observationally preferred spectral index n = 1 
(/3 = -2,7 = 0). 



VII. WHAT INFLATIONARY THEORY SAYS ABOUT DENSITY PERTURBATIONS, 
AND WHAT SHOULD BE SAID ABOUT INFLATIONARY THEORY 

For many years, inflationists keep insisting that in contrast to the generation of gravitational 



waves 



3, 



which begins with the amplitude h ~ Ipi/X ("half of the quantum in the mode") and 
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finishes with the amplitude h ~ lpi/{c/H) ~ H/Hpi and power /i^ ~ (H/Hpi)'^ (see ^ and Sees. 
mil lIVp . the generation of primordial density perturbations (i.e. scalar metric perturbations) is 
more efficient by many orders of magnitude. Contrary to the calculations reviewed in Sec. IVH 
inflationary theory claims that in a cosmological model with the same H, the resulting amplitude 
and power of the quantity (" should contain a huge extra factor, tending to infinity in the limit of 
the standard de Sitter inflation. 

Traditionally, the claimed inflationary derivation of density perturbations goes along the fol- 
lowing lines. One starts from the ground state quantum fluctuations taken from a different theory, 
namely from a theory of a free test scalar field, where metric perturbations are ignored altogether, 
and writes 

Then, from the gauge transformation (jSip . assuming that the l.h.s. is equal to zero, one finds the 
characteristic time interval 5t "to the end of inflation" and puts into 5t the estimate of 5ip from 
the above-mentioned 'quantum' evaluation: 

Then, one declares that the dimensionless ratio 6t/H^^ is what determines the dimensionless metric 
and density variation amplitudes in the post-inflationary universe. And this is being presented as 
the "famous result" of inflationary theory: 

p 2-nipQ ipo/H' 

This "inflationary mechanism" of generation of density perturbations is claimed to have been 
confirmed in numerous papers and books 



58| . The most dramatic feature of the claimed inflationary 
result is the factor <pq/H in the denominator of the final expression. I call it a 'zero in the 
denominator' factor [15]. According to Eq. (14ip this factor is ^77; in the literature, this factor 
appears in several equivalent incarnations, including such combinations as Vip/V, 1 +p/ pc^ = 
1 + w, e, where e = 7 = —H/H^, and so on. 

The early Universe Hubble parameter H featuring in the numerator of the final expression 
cannot be too small because H{t) in scalar field driven cosmologies is a decreasing, or at most 
constant, function of time. The parameter H in the very early Universe cannot be smaller than, 
say, the H in the era of primordial nucleosynthesis. So, the numerator of the "famous result" 
cannot be zero, but the denominator can, at 7 = 0. Therefore, the "famous result" prescribes the 
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arbitrarily large numerical values to the amplitudes of density perturbations in the limit of the 
de Sitter evolution 7 = (/? = —2). One has to be reminded that it is this gravitational pump 
field that generates gravitational waves and other perturbations with primordial spectrum n = 1, 
advocated long ago on theoretical grounds by Harrison, Zeldovich, Peebles, and which is in the 
vicinity of the primordial spectral shape currently believed to be preferred observationally. 

The inflation theory claims to have predicted a "nearly" scale-invariant spectrum of density 
perturbations because on the strictly scale-invariant spectrum the predicted amplitudes blow up 
to infinity. 

More recent literature operates with the 'curvature perturbation TZ\ equivalent to our (" from 
Eq. (|T7|) . A typical quotation states: "The amplitude of the resulting scalar curvature perturbation 
is given by 

where H is the Hubble parameter, cp is the time derivative of the inflaton cj), and S(p is the infla- 
ton fluctuation on a spatially flat hypersurface. The quantum expectation value of the inflaton 
fluctuations on super-horizon scales in the de Sitter space-time is 




Here again one is invited to believe that the ground state quantum fluctuations in one theory 
(without metric perturbations) are responsible for the appearance of arbitrarily large (factor 7 
in the denominator of the TZ power spectrum) gravitational field fluctuations in another theory. 
Inflationists are keen to make statements about the resulting gauge-invariant curvature perturba- 
tions without having included metric perturbations in the initial conditions. The infinitely large 
curvature perturbation TZ is supposed to occur at 7 = 0, that is, exactly in the de Sitter model 
from which all the 'quantum' reasoning about scalar field fluctuations has started. 

According to inflationary views on the generating process, the amount of the created scalar 
particles-perturbations is regulated not by the strength of the external gravitational 'pump' field 
(basically, space-time curvature in the early Universe) but by the closeness of the metric to a de 
Sitter one. In a sequence of space-times with very modest and approximately equal values of H 
you are supposed to be capable of generating arbitrarily large amplitudes of the scalar metric 
perturbation by simply going to smaller and smaller values of H. 

Since the quantity TZ (as well as the quantity C and the gravitational wave function h) oscillates 
with a slowly decreasing amplitude in the initial short-wavelength regime (see Eq. ()48p ) and remains 
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constant during the subsequent long-wavelength regime, the arbitrarily large amplitude of the 
resulting "inflation-predicted" scalar curvature perturbation TZ must have been implanted from 
the very beginning, i.e. from the times in the high-frequency regime of evolution, when the intial 
quantum state for density perturbations was defined. It is important to remember that for many 
years inflationists claimed that the reason for the huge difference between the resulting scalar 
and tensor perurbations was the "big amplification during reheating" experienced by the long- 
wavelength scalar metric perturbations. These days, the explanation via the "big amplification 
during reheating" is not even mentioned. These days, the most sophisticated inflationary texts put 
forward, as the foundation for their belief in arbitrarily large resulting scalar metric perturbations 
TZ and the "Bunch-Davies vacuum", i.e. a concept from the theory of a test scalar field, where 
the gravitational field (metric) perturbations are absent altogether. (Moreover, the Bunch-Davies 
vacuum was originally introduced as a de Sitter invariant state prohibiting the particle production 
by definition.) The absurd proposition of inflationists with regard to the density perturbations 
is sometimes called a "classic result", and it is widely used for derivation of further incorrect 
conclusions in theory and data analysis. 

The most common, and most damaging, inflationary claim is that the amplitudes of relic grav- 
itational waves must be "suppressed", "sub-dominant", "negligibly small" in comparison with 
primordial density perturbations, especially for models with 7 — )• 0. Having arrived at a divergent 
formula for density perturbations, i.e. with a 'zero in the denominator' factor, inflationists compose 
the ratio of the gravitational- wave power spectrum Vt to the derived divergent scalar metric power 
spectrum (the so-called 'tensor-to-scalar ratio' r = Vt/'P(^) and write it as 



The inflationary theory binds the amplitude of scalar perturbations with the spectral index, and 
makes the absurd prediction of arbitrarily large amplitudes of density perturbations in the limit 
of models with e = 7 = (n^ = 1, nt = 0). But the inflationary "consistency relation" (j66p 
encourages and misleads one to believe that everything is perfect - as if it were the amount of 
gravitational waves that must go to zero in this limit. To make the wrong theory look acceptable, 
inflationary model builders keep "P^ fixed at the observationally required level of 10~^ or so, and 
move H/Hpi down whenever e goes to zero in the inflationary divergent formula 



r = 16e = —Snt- 



(66) 




thus making the amount of relic gravitational waves arbitrarily small. 
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It is instructive to see how the problem of initial conditions is delt with by S. Weinberg [27]]. The 
author operates with equations for ( (equivalent to our Eq. (^5]) ). gravitational waves (equivalent 
to our Eq.dJ]) for h), and a test scalar field a (which is known [2] to be similar to the equation 
for gravitational waves ). When it comes to the initial conditions, the author says that they are 
"designed to make" (^ipo/H, h, and a behave like conventionally normalized free fields in the remote 
past. In other words, the initial conditions for h and a are not "designed to make" a potentially 
vanishing factor to enter the normalization, but the initial conditions for (" are. The normalization 
of ( becomes now propotional to l/\/e, and the power spectrum of ( acquires the factor e in 
the denominator. The author calls this divergent power spectrum of scalar metric perturbations a 
"classic" result. 

One can imagine why this divergent formula is called "classic". Something repeated in the 
literature so many times could become "classic" more or less automatically, regardless of its true 
value. However, Weinberg has not explicitly stated that the "classic" result is a correct result. 
On the contrary, the recent paper 2g], which deals with gravitational waves, makes assumptions 



diametrically opposite to the prescriptions of the "classic" result (but the paper does not explicitly 
say that the "classic" result is an incorrect result). That paper chooses for the analysis a model 
with rit = and 'tensor /scalar ratio' r = 1. This choice is in agreement with conclusions of 
the quantum theory that I advocated and reviewed in Sec lVH but it is in confiict with what is 
demanded by infiationary Eq. (l66p (quoted also in [29|). Indeed, the "classic" result for the case 
Tit = implies the non-existence of the very subject of discussion, namely, relic gravitational waves. 
Although the choice = 0, r = 1 has been made [28 ] purely for numerical convenience, it seems to 



me that the cautious formulation "designed to make" 



recent 



27l |, together with the earlier 29| and more 



28t | treatments, testify to a certain evolution of views on the subject of initial conditions. 



It appears that some attempts of technical derivation of the infiationary (^-normalization suffer 
from serious inaccuracies in dealing with quantum operators and quantum states. Certainly, it is 
incorrect to think that by "demanding that and a obey the standard creation and annihilation 
commutation relations we get a normalization condition for ^". To put it in the context of a 
medical analogy: the rules for a surgical operation do not identify the patient on whom you want 
to operate. Let us discuss this point in more detail. 

Suppose, being (mis)guided by various infiationary prejudices, you decided to write, instead of 
Eqs. (l59p . (I60p . the following asymptotic expressions for the operators q, p: 



hao 1 
2T7f 



j^g-m(7?-7?o) _|_ l3tg«"{»?-'?o) 



(67) 
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p = ij -—^ r_be-*"(^-^o) + bte^''(''-''o)l . (68) 
V 2 ao L J 

The commutation relationships [q, p] = ih dictate exactly the same commutation relationships for 
the operators b^, b as they did for the operators c^, c. Namely, [b, b^^] = 1. 

The commutation relationships for the annihilation and creation operators are exactly the same, 
but the quantum state |0s) annihilated by b, 

b|0,) = 0, 

is totally different, it is not the ground state of the Hamiltonian (j54p . Calculation of the mean 
square value of the variable ^ and its canonically-conjugate momentum gives at rj = r]Q: 

(0,|q>,) = (0,|p2|0,) = ^70, 

2 70 2 

and the factor ^Jjq cancels out in the uncertainty relation: 

AqAp = ^. 

The initial rms value of C is proportional to I/^/jq and therefore contains the 'zero in the 
denominator' factor, but this happens only because the quantum state |0s) is an excited (multi- 



quantum) squeezed vacuum state (for more details, see [15| and references there). The choice of 
this state as an initial state for (^-perturbations would make them arbitrarily large, in the limit 
of 7o — )• 0, right from the very beginning, i.e. from the time of imposing the initial conditions at 
r] = r]Q. The multi-quantum state \0s) is not a choice of the initial state that is regarded physically 
motivated. 

At the level of classical equations, inflationists have reacted to the visual analogy between 
equations (JI]), (f38l) . rather than, say, to the visual analogy between equations ([2]), (H8]l . They 
have modeled initial conditions for the function fis on the initial conditions for the function /i-p, 
instead of modeling initial conditions for on the initial conditions for the g.w. function h. The 
inflationary initial conditions are usually written in the form 

^ —ikr) 



1,^T,S) = -^^^ ' for k\r]\ » 1. 



These initial conditions are correct for gravitational waves, but are incorrect for density perturba- 
tions. As is seen from Eq. ()47p . these initial conditions for fis would require the gauge-invariant 
metric perturbation ^ (as well as curvature perturbation TZ) to be arbitrarily large, in the limit of 
models with 7 — )■ 0, right from the very beginning, i.e. from the time of imposing initial conditions 

33 



in the short-wavelength regime. In other words, the inflationary initial conditions are "designed to 
make" the quantity C to be divergent as 1/^/7 from the start. These incorrect initial conditions 
are used in all inflationary calculations and conclusions, including the latest claims of inflationists 
on what the observations "really" tell them about inflation. 

Finally, let us assume that the 'correct' Lagrangian is given by Eq. (l57p . (Surely, you may 
assume correct whichever Lagrangian you wish, because it is taxpayers [30] who will be paying 
the price at the end of the day.) In contrast to the Lagrangians (jlSp . (j55p . this new Lagrangian 
vanishes in the most interesting limit of models with 7—7-0. 

For this new Lagrangian, the canonical quantization would require us to write (instead of 
Eqs.dSl]), dSOD): 

'/iOoVTo 1 



q 



2 0^/7 \/7o 
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(69) 



(70) 



where [q, p] = ih and [d, d^^] = 1. 

The ground state \^new) of the new Hamiltonian associated with this new Lagrangian obeys the 
condition 

d|0„e«,) = 0. 

The mean square value of the variable in this state at 7/ = 770 is given by 

I 2| 

2 70 

Technically speaking, one could argue that although the initial rras value of C is divergent as IZ-^tq, 
the divergency takes place for the ground state (of this new Hamiltonian), not for an excited state 
(of this new Hamiltonian). But this technical subtlety requires a vanishing Lagrangian, and in any 
case it does not change much from the physical point of view. 

Indeed, independently of technical arguments, it is important to realize that a proposal for a 
divergent scalar metric power spectrum in the limit of n = 1, whatever the reasons for this proposal 
might be, is in conflict with available observations. The currently derived 'best fit' value for n is 
slightly lower than 1, but the data allow n = 1, even if with a smaller likelihood. This means 
that for the tested spectral indeces in the vicinity and including n = 1 the data are consistent 
with finite and small amplitudes at the level of the best fit amplitude (as implied by the quantum 
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theory discussed in Sec lVip . But there is absolutely nothing in the data that would suggest the need 
for a catastrophic growth of the amplitude (demanded by inflationary theory) when the data are 
fitted against spectra with indeces approaching and crossing n = 1. (Surely, the same comparison of 
inflationary predictions with available observations is addressed by inflationists as "one of the most 
remarkable successes ... confirmed by observations".) What would we conclude about the assumed 
Lagrangian (j57p whose associated Hamiltonian leads to predictions contradicting observations ? 



We "would conclude that it was the wrong Lagrangian" 

The inflationary 'zero in the denominator' factor (if you decided to commit suicide and include 
it in the scalar metric power spectrum) should be taken at the moment of time, for a given mode n, 
when the mode begins its superadiabatic evolution. In general, this factor is n-dependent, and hence 
it affects not only the overall normalization of the spectrum, but also its shape. It enables one to 
'generate' a flat, or even a blue, power spectrum of scalar perturbations by gravitational pump fields 
which in reality can never do this. It enables one to derive all sorts of wrong conclusions in various 
subjects of study, ranging from the formation of primordial black holes and up to perturbations in 
cyclic and brane-world cosmologies. In particular, recent claims stating that a given ultra-modern 
theory predicts an "extremely small r < 10"^'* " or such an r that "a tensor component. ..is far below 
the detection limit of any future experiment", mean only that the predictions were based on the 
incorrect (inflationary) formula for scalar metric perturbations, with the 'zero in the denominator' 
factor. Whatever the incorrectly derived r may be, small or large, the use of inflationary 



in data analysis (what, unfortunately, is regularly being done in the CMB data analysis 
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33]) can only spoil the extraction of physical information on relic gravitational waves from the 



data. 

The self-contradictory nature of conclusions based on the inflationary formula (j66p is dramat- 
ically illustrated by the claimed derivation of limits on the amount of gravitational waves from 



the WMAP data 



34l |. The authors use Eq. (j66p and explicitely quote, in the form of their equa- 



tions (17) and (18), the power-spectrum amplitude of inflationary scalar perturbations and the 
'tensor-to-scalar' ratio r: 

- ^ r - 166 fl) 

The conclusion of this (and many other subsequent papers) is such that according to the likelihood 
analysis of the data the value of r must be small or zero. In other words, the maximum of the 
likelihood function is either at r = or, in any case, the value of the observed r is perfectly well 
"consistent with zero". Comparing this conclusion with the inflationary formulas (I), quoted in 
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the same paper, one has to decide whether this conclusion means that the WMAP data are also 
perfectly well consistent with arbitrarily large scalar amplitudes A^, and hence with arbitrarily 
large CMB anisotropics caused by density perurbations, or that the cited inflationary formulas 
are wrong and have been rejected by observations (see also a discussion below, in Sec lXIip . The 
CMB data and the inflationary theory of density perturbations are in deep conflict with each other 
for long time, but inflationists and their followers keep claiming that they are in "almost perfect 
agreement" . 

It seems to me that the situation in this area of physics and cosmology remains unhealthy for 
more than 25+ years. (This is my mini-version of 'an obligation to inform the public' - from 



35l|.) It appears that more than 2+ generations of researchers have been 'successfully' misled 



by the "standard inflationary results". One can only hope that the present generation of young 
researchers will be smarter and more insightful. 



VIII. WHY RELIC GRAVITATIONAL WAVES SHOULD BE DETECTABLE 

The detailed analysis in previous sections is crucial for proper understanding of the very status 
of relic gravitational waves. Are we undertaking difficult investigations because we want to find 
an optional "bonus" , "smoking gun" , "limitation" on dubious theories, or we search for something 
fundamental that 'must be there', and at a measurable level ? 

Strictly speaking, it is still possible that the observed CMB anisotropies have nothing to do 
with cosmological perturbations of quantum-mechanical origin, that is, with the superadiabatic 
evolution of the ground state of quantized perturbations. The first worry is that, even if the su- 
peradiabatic (parametric) mechanism is correct by itself, the 'engine' that drove the cosmological 
scale factor was unfortunate and the pump field was too weak. In this case, a relic gravitational- 
wave signal (as well as scalar perturbations of quantum-mechanical origin) could be too small for 
discovery. I think this possibility is unlikely. First, it is difficult to imagine an equally unavoidable 
mechanism - basic laws of quantum mechanics and general relativity - for the generation of the 
presently existing, as well as the processed in the past, long-wavelength cosmological perturba- 
tions. Second, the inevitable quantum-mechanical squeezing and the standing-wave character of 
the generated perturbations, which, among other things, should have resulted in the scalar metric 
power spectrum oscillations (see Sec lIVp . seems to have already revealed itself {s^] in the observed 
CMB power spectrum oscillations. (These oscillations are often being associated - in my opinion, 
incorrectly - with baryonic acoustic waves and are called "acoustic" peaks). 
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A second worry is that we may be wrong in extrapolating the laws of general relativity and 
quantum mechanics to extreme conditions of the very early Universe. Although we apply these 
laws in environments that are still far away from any Planckian or 'trans-Planckian' ambiguities, 
it is nevertheless an extremely early and unfamiliar Universe (which we want to explore). In 
principle, it is possible that something has intervened and invalidated the equations and rules that 
we have used for derivation of relic gravitational waves. This would also invalidate the equations 
and rules used in the derivation of density perturbations, but the generation of scalar perturbations 
requires an extra hypothesis in any case. In other words, even if the driving cosmological 'engine' 
was right, the employed quantum theory of arising perturbations could be wrong. Hopefully, this 
complication also did not take place. 

Assuming that the dangers did not materialize, we come to the conclusion that we have done 
our job properly. The theoretical calculations were adequate, the normalization of hrms on the 
CMB lowest-order multipoles was justified, and therefore the relic gravitational wave background 
must exist, and probably at the level somewhere between the dotted and solid lines in Figj5j The 
other side of the same argument is that if we do not detect relic gravitational waves at this level, 
something really nasty, like the above-mentioned dangers, had indeed happened. From not seeing 
relic gravitational waves, we would at least learn something striking about the limits of applicability 
of our main theories. 

It is important to stress again that what is called here the relic gravitational waves is not the 
same thing which is sometimes called the inflationary gravitational waves. As the name suggests, 
statements about the inflationary gravitational waves are based on the inflation theory, as applied to 
density perturbations and gravitational waves. Although, conceptually, the inflationary derivation 
of density perturbations is an attempt of using the mechanism of superadiabatic (parametric) 
amplification (see SecH]), the actual implementation of this approach has led inflationists to the 
divergency in the scalar metric spectrum (see Sec iVIip . which has been converted, for the purpose 
of "consistency", into a statement about small r. In particular, for the "standard inflation" (7 = 
e = 0), the inflationary theory predicts (see Eq. (|66p ) the "smoking gun of inflation" in the form 
of a zero amount of inflationary gravitational waves (r = 0). This is in sharp contrast with the 
superadibatic (parametric) prediction of a finite and considerable amount of relic gravitational 
waves, as shown by a dotted line in Figl5j 

Some recent papers, still based on the incorrect theory with Eq. (l66p . start making additional 
artificial assumptions about the "natural" inflationary conditions, which amounts to postulating 
that e should not be too small, but should, instead, be at the level of, say, 0.02. This kind of 
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argument gives some authors "more reason for optimism" with regard to the detection of primordial 
gravitational waves. Surely, an incorrect theory can lead to predictions not very much different 
from predictions of a correct theory in some narrow range of parameters that are supposed to 
be deduced from observations. One should remember, however, that theories are tested not only 
by what they predict but also by what they do not predict. Any observation consistent with the 
parameter value other than you postulated is against your theory. 

Since the spectrum of relic gravitational-wave amplitudes is decreasing towards the higher fre- 
quencies, see Figini it is the lowest frequencies, or better to say the longest wavelenghts, that 
provide the most of opportunities for the (indirect) detection. It is known for long time 
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that gravitational waves affect the CMB temperature and polarization. The low-order CMB mul- 
tipoles {i < 100) are mostly induced by cosmological perturbations with wavelengths ranging from 
10 times longer and up to 10 times shorter than the present-day Hubble radius Ih- And this is the 
range of scales that will be in the center of our further analysis. 



IX. INTENSITY AND POLARIZATION OF THE CMB RADIATION 

A radiation field, in our case CMB, is usually characterized by four Stokes parameters 
(/, Q,U,V) [l^. / is the total intensity of radiation (or its temperature T), Q and U describe 
the magnitude and direction of linear polarization, and V is the circular polarization. The Stokes 
parameters of the radiation field arriving from a particular direction in the sky are functions of a 
point 6,(f) on a unit sphere centered on the observer. The metric tensor Qat on the sphere can be 
written as: 

da^ = Qabdx^dx^ = dO'^ + sin^ Odct)'^. (71) 

The Stokes parameters are also functions of the frequency of radiation but angular dependence 
is more important for our present discussion. 

The Stokes parameters are components of the polarization tensor Pah [3] which can be written 

I ( I + Q - (U-iV)sm9 \ 

^\-{U + iV)sme {I-Q)sm^eJ 

As every tensor, the polarization tensor Pat transforms under arbitrary coordinate transformations 
on the sphere, but some quantities remain invariant. 

We can build invariants, linear in Pab and its derivatives, with the help of the tensor gab and 
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the unit antisymmetric pseudo-tensor e°^. First two invariants are easy to build: 



1(9, 



(73) 



Two other invariants involve second covariant derivatives of Pab [ITl] : 



-2{Pab{e,4>)r\ B 



-2{Pab{9A)y'"'e'^d, 



(74) 



Being invariant, quantities B do not mix with each other. To calculate E and i? in a given 
point 0, we do not need to know the polarization pattern all over the sky, but we do need to know 
the derivatives of Pab at that point. Whatever the numerical values of or i? in a given point are, 
E and B will retain these values under arbitrary coordinate transformations (smoothly reducable 
to an ordinary rotation). The invariant B is built with the help of a pseudo-tensor and 
therefore i? is a pseudo-scalar rather than an ordinary scalar. While E does not change sign under 
a coordianate reflection, B does. With B one can associate the notion of chirality, or handedness 
(compare with polarization tensors (jl2p . (jlSp ). Clearly, if given cosmological perturbations are 
such that they themselves are incapable of supporting the handedness, it will not arise in the CMB 
polarization which these perturbations are responsible for. 

The invariant quantities /, V , E, B can be expanded over ordinary spherical harmonics 



v{eA) 
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(75a) 
(75b) 
(75c) 
(75d) 



The set of (complex) multipole coefficients {aj^, a^^^, af^, af^) completely characterizes the inten- 
sity and polarization of the CMB. 

The same multipole coefficients participate in the expansion of the tensor Pab itself, not only in 
the expansion of its invariants, but the expansion of Pab requires the use of g eneralized spherical 
functions, the so-called spin- weighted or tensor spherical harmonics 42, [43] • For example, the 
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tensor Pat can be written as 

^ oo £ 

£=0 m=-e 
^ oo I 



where 



and 5^(L)a6( 



=2 m=-« 

are the 'gradient' and 'curl' tensor spherical harmonics 



In what follows, we will not be considering the circular polarization V ^ and we will sometimes 
denote the multipole coefficients collectively by a^, where X = I,E,B. These coefficients are to 
be found from solutions to the radiative transfer edquations in a slightly perturbed universe. 



X. RADIATIVE TRANSFER IN A PERTURBED UNIVERSE 

Polarization of CMB arises as a result of Thompson scattering of the initial 
on free electrons residing in a slightly perturbed universe, Eq.dS]). Following 



y unpo. 



4l|, 



arized light 



4j], [40|, it is 



convenient to describe Stokes parameters in terms of a 3-component symbolic vector n: 



n 



\n3j 



I-Q 
-2U 



In the zero-order approximation, we assume that all hij 
fully homogeneous, isotropic, and unpolarized. Then, 

n'-'^^ = no(z^a(r/))u. 



(76) 



and the CMB radiation field is 



(77) 



where 



u 



1 



In the presence of metric perturbations, we write 



(78) 



where n^^^ is the first order correction. The functions n*-^-* depend on {rj, x*, u, e*), where u = ualrj) 
and e* is a unit spatial vector along the photon's path. Our final goal is to predict, with as much 
completeness as possible, the values of the Stokes parameters at the time of observation r] = r]R. 
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The functions n^^^ satisfy the hnear version of the radiation transfer equations. These equations 
can be written [ij 



f{i>)no{i>) i jdhij 

e e-' — - — 

2 drj 
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(79) 



The astrophysical inputs from unpolarized radiation and free electrons are described, respectively, 
by /(P)no(t') and q{r]), while the scattering process is described by the Chandrasekhar matrix 
P(e*; e'^). The input from the gravitational field (metric) perturbations is given by the combination 
dhij / drj . Certainly, when all hij = 0, all n^^^ vanish if they were not present initially, which we 
always assume. 

The combination e^e^dhij/dr] gives rise to disparate angular structures for gravitational waves 
and density perturbations. Let us consider a particular Fourier mode with the wavevector n = 
(0, 0, n). The polarization tensors of gravitational waves generate the structure [l^] 

eV4(n) = (l-//2)e±2^<A^ (80) 

where /i = cos and the it signs correspond to the left s = 1 = L and right s = 2 = R polarization 
states, respectively. Solving Eq. (j79p for n'-''"^ in terms of a series over e*™"^, one finds that the terms 
proportional to e^^**^ are retained while other components e*™''^ are not arising. This is because 
the Chandrasekhar matrix does not create any new m(p dependence. 

In contrast to gravitational waves, the same metric combination with polarization tensors (jl3p 
of density perturbations produces the structure which is (/>- independent. Although in the case of 
density perturbations, Eqs. ()79p contain one extra term, proportinal to the electron fluid velocity, 
e*fi = —ifiVb (expressible in terms of metric perturbations via perturbed Einstein equations), this 
term is also (/>- independent. Since the invariant B depends on the derivative of n^-*") over (p, one 
arrives at the conclusion that B = for density perturbations and B ^ for gravitational waves. 
These results for one special Fourier mode n = (0, 0, n) can then be generalized to any arbitrary 
n. As we anticipated on the grounds of handedness, gravitational waves can generate the i?-mode 
of CMB polarization, but density perurbations can not 
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XI. STATISTICS AND ANGULAR CORRELATION FUNCTIONS 



The linear character of the radiation transfer equations (|79p makes the randomness of hij being 

inhereted by n^^^. A consistent handhng of Eq. (|79p requires us to use the spatial Fourier expansion 

+00 

-00 s=l,2 

(81) 

where random coefficients are the same entities which enter the metric field Eq.([9]). The CMB 
anisotropies are random because the underlying metric perturbations are random. 

In the strict quantum-mechanical version of the theory, where c^. are quantum-mechanical op- 
erators, the CMB field n^^^ itself becomes a quantum-mechanical operator. If the initial quantum 
state of the system is chosen to be the ground state |0), Cn |0) = 0, all statistical properties of the 
system are determined by this choice. 

Note the extra degree of uncertainty that we will have to deal with. Usually, the signal is 
deterministic, even if totally unknown, and the randomness of the observed outcomes arises at the 
level of the measurement process, as a consequence of the uncontrollable noises. In our problem, if 
cosmological perturbations do indeed have quantum-mechanical origin, the signal itself is inherently 
random and is characterized by a quantum state, or a wave-function, or a probability distribution 
function. At best, we can predict only a probability distribution function for possible CMB maps 
(outcomes). This is true even if the dynamics and cosmological parameters are strictly fixed and 
the measurement process is strictly noiseless. 

Each mode of cosmological perturbations has started its life in the initial vacuum state. This 
state is Gaussian. In course of time it evolved into a squeezed vacuum state. Squeezed vacuum 
states retain the Gaussianity, even though developing the strongly unequal variances in amplitudes 
and phases. Therefore, the actually observed coefficients (our own realisation of the CMB map 
belonging to the theoretical ensemble of CMB maps) are supposed to be drawn from the zero-mean 
Gaussian distributions for a^. If the observed CMB map looks a bit strange to you, due to the 
presence, for example, of some hints on 'axes' and 'voids', this is not necessarily an indication of 
non-Gaussianity of the underlying cosmological perturbations. Even if the observed CMB map 
consists entirely of your own images, this is not a proof of non-Gaussianity. It is a legitimate 
procedure to postulate some sort of a non-Gaussian distribution, by introducing a new parameter 
/tvl, and try to find fiy^ from the data of a single actually observed map. But it is even more 
a legitimate procedure to insist on Gaussianity of perturbations, because it follows from the very 
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deep foundations explained above, and therefore regard f nl = hy definition. From this position, 
if the set of the observed coefficients af^ looks a bit strange to you, this is simply because our 
Universe is not as dull and 'typical' as you might expect. 

In this paper, we simplify the problem and treat cosmological fields classically rather than 
quantum-mechanically. We also make mild statistical assumptions. Specifically, we assume that 
classical random comples numbers Cn satisfy a limited set of statistical requirements: 

(Cn) = {Cn') = 0, {Lcn') = (Ccn') = Sss'S^^\n - n'), {Lin') = {cjn') = 0, (82) 

where the averaging is performed over some probability distributions. We do not even assume 
outright that these distributions are Gaussian. The rules (j82p are sufficient for the most of our 
further calculations. 

We want to know the value of quantities at the time of observation r] = rjji. To find these 
quantities we have to integrate Eq. (j79p over time, with all the astrophysical and gravitational inputs 
taken into account. The derived coefficients af^ are random, because the participating coefficients 

are random. We can calculate various correlation functions of the CMB by calculating the 
quantities (a^a^^,). Using the rules ([82]) . one can show that these quantities take the form 

{af^*a^^,) = Cf^ Su'Smm', (83) 

where Cf^ depend on the gravitational mode functions and astrophysical input. Cf'^ are cal- 
culable as general theoretical expressions. 

The quantities Cf^' are called the multipoles of the corresponding CMB power spectrum XX' . 
As usual, the power spectrum of a field contains less information than the field itself, but we will 
mostly ignore this loss of information. Also, it is worth noting that practically any feature in 
the CMB power spectrum can be "predicted" and "explained" solely by the properly adjusted 
primordial spectrum of cosmological perturbations. But we will not go along this line and will 
stick to simple power-law primordial spectra. 

In the case of relic gravitational waves, the general shape of today's CMB power spectra and their 
features in the £-space are almost in one-to-one correspondence with the processed metric power 
spectra and their features in the wavenumber n-space. The processed metric power spectra should 
be taken at the time of decoupling of CMB. The TT power spectrum is determined by the power 
spectrum of the metric itself, hh. The EE and BB power spectra are largely determined by the 
power spectrum of the metric's first time-derivative, h'h' . And TE power spectrum is determined 
by the cross power spectrum of the metric and its first time-derivative, hh' . For the gravitational- 
wave background with parameters indicated by dotted line in Figj5l the corresponding metric and 
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FIG. 7: The left panel shows (a) the power spectrum of temperature anisotropics £{i + l)Cj^ (in /iK^) 
generated by (b) the power spectrum of g.w. metric perturbations hh (|29p . /? = —2. The right panel shows 
(c) the power spectra of polarization anisotropics l{l + l)Cf^ (solid line) and £{£ + l)Cf^ (dashed line), 
panel (d) shows the power spectrum of the first time derivative of the same g.w. field, h'h' . 

CMB power spectra are shown in FiglT] and FiglSj For the gravitational-wave background with 
parameters indicated by a solid line in FiglU the summary of CMB spectra is shown in FiglH The 
summary also includes the reionization 'bump' at i < 12. (More details about these figures can be 
found in [17]). 

All CMB power spectra have been calculated as averages over a theoretical ensemble of all 
possible realizations of the CMB field. In its turn, this randomness of CMB anisotropics ensues 
from the randomness of the gravitational field coefficients c^, c^- The characteristic parameters 
of a stochastic field, such as its mean values and variances, are, by definition, averages over the 
ensemble of realizations. However, in CMB observations, we have access to only one realization of 
this ensemble, which can be thought of as a single observed set of coefficients a^. Is it possible to 
find the parameters of a stochastic process by studying only one realization of this process ? 

The answer is yes, if the correlations of the stochastic process decay sufficiently quickly at large 
separations in time or space where the process is defined [47!]. The process allowing the derivation 
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FIG. 8: The bottom panel shows the cross-po'wer spectrum hh' of gravitational "waves, "whereas the top 
panel sho"ws the angular po"wer spectrum £{£ + l)Cj^ caused by these -waves. The negative values of these 
functions are depicted by broken lines. 

of its true parameters, "with probability arbitrarily close to 1, from a single realization is called 
ergodic. For example, the distribution of galaxies, or a stochastic density field, in an infinite 3 - 
space (our Universe) may be ergodic, and then by studying a single realization of this distribution 
we could extract the true parameters of the underlying random process. In the theory that we 
are discussing here, the randomness of the linear density field is also described by the random 

S S * 

coefficients Cn, appearing in the metric perturbations. 

If the process is non-ergodic, there "will be an inevitable uncertainty surrounding the parameter's 
estimation derived from a single realization. This uncertainty should not be mis-taken for the 
'cosmic variance', often quoted and plotted on the observational graphs. The cosmic variance is a 
mathematically correct statement about the size of the variance of the x^-distribution for 2i + 1 
independent Gaussian variables. There is nothing particularly cosmic in the cosmic variance. 
The elementary theorem, called cosmic variance, has nothing to say about the ergodicity or non- 
ergodicity of a given stochastic process. In particular, the universal validity of cosmic variance 
cannot prevent the extraction of exact parameters from the observation of a single realization of 
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FIG. 9: The summary of CMB temperature and polarization anisotropies due to relic gravitational waves 
with n = 1.2 and R = 1. 

the stochastic density field in our Universe, if the density field is ergodic, - and it might be ergodic. 

The problem is, however, that on a compact 2-sphere, where the random CMB is defined, 
ergodic processes do not exist. We will always be facing some uncertaintly related to non-ergodicity. 
The size of this uncertainty about the derived parameter depends on the statistics and employed 
estimator. Under some conditions, this uncertainty is close, numerically, to the size of the usually 
quoted cosmic variance [4^. This discussion is important, because we are now approching the 



observational predictions and the ways of discovering relic gravitational waves in the CMB data. 



XII. TEMPERATURE-POLARIZATION CROSS-CORRELATION FUNCTION 

The numerical levels of primordial spectra for tensor and scalar metric perturbations are approx- 
imately equal. This means that the amplitudes of those Fourier modes n which have not started 
yet their short-wavelength evolution are numerically comparable, and the metric mode functions 
are practically constant in time. In contrast, in the short-wavelength regime, the amplitudes of 
gravitational waves adiabatically decrease, while the amplitudes of gravitational field perturbations 
associated with density parturbations may grow. Specifically, at the matter-dominated stage, the 
function hi(rj) grows and overtakes h(rj), which remains constant. 

In the context of CMB anisotropies, the crucial time is the epoch of decoupling. The wavelenghts 
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FIG. 10: The dotted line shows the contribution of density perturbations alone, while the dashed line shows 
the contribution of gravitational waves alone. The solid line is the sum of these contributions. It is seen 
from the graph that the inclusion of g.w. makes the total curve to be below the d.p. curve. 



of modes with n < 100 were comfortably longer than the Hubble radius at the decoupling. The 
influence of these modes, both in gravitational waves (g.w.) and density perturbations (d.p.), have 
been projected into today's XX' anisotropics at £ < 100. In this interval of i, the g.w. contribution 
is not a small effect in comparison with the d.p. contribution. The g.w. contribution to the CMB 
power spectra is illustrated in FiglTl FiglH and FigiH 

One way of detecting relic gravitational waves is based on measuring the BB auto-correlation. 
This method is clean, in the sense that density perturbations do not intervene, but the expected 
signal is very weak, see FiglJl FiglD) 

We propose [17| to concentrate on the TE cross-correlation (without, of course, neglecting the 



BB searches). The TE signal is about two orders of magnitude stronger than the BB signal, and 
the use of a cross-correlation is always better than an auto-correlation, in the sense of fighting the 
noises. The special feature allowing to distinguish the g.w. part of TE from the d.p. part of TE 
is the difference in their signs. In the interval i < 100 the {TE)gw must be negative (see FigJS]), 
while the {TE)dp must be positive at lowest £'s and up to, at least, £ ^ 50. This difference in 
sign of TE correlation functions is the consequence of the difference in sign of the g.w. and d.p. 
metric cross-power spectra hh' in the interval n < 70 \v\. (The difference in sign of QT correlation 
functions is discussed in the earlier paper [49].) 

An example of expected g.w. and d.p. contributions to the TE correlation function is shown 
in FigllOi To include reionization and enhance the lowest i's, we plot the function {e + l)/2TT Cf^ 
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rather than the usual £{i + l)/27r Cj ■ The g.w. and d.p. metric power spectra are normahzed in 
such a way that they give R = 1, where 

More specificahy, they give equal contributions to the total 'best fit' [32] temperature quadrupole 
C^^. The dashed line in FigJlOl shows the effect of the g.w. background marked by dotted line 
in Figl21 The dotted line in FigllOl shows the effect of d.p. metric perturbations, with the same 
primordial spectral index, n = 1. The sum of the two contributions is shown by a solid line. The 
fact that the two contributions may almost cancel each other does not mean that the g.w. signature 
is weak and hard to measure. The g.w. signal is in the strong deviation of the total TE spectrum 
from the expectation of the d.p. model. 

It is important to remember that the widely publicized "tight limits" on the 'tensor-to-scalar 
ratio' r which allegedly rule out R = 1 and any R ^ 0, unless R is small, were derived by 
making use of the inflationary "consistency relation" (see Sec lVIip . In these derivations, the g.w. 
spectral index rit is being taken from the relationship nt = —r/8, which automatically sends r to 
zero when nt approaches zero. If this relationship is regarded as an artificial extra condition on 
g.w. parameters, then the results of such a data analysis may be of some value to those who are 
interested in this ad hoc condition, but not to those who are interested in determination of the true 
amount of relic gravitational waves. If, on the other hand, this relationship is regarded as part of 
inflationary theory, then such a data analysis is deeply self-contradictory. Indeed, the invariably 
derived conclusion, according to which the maximum of the likelihood function for r is at r = 0, 
or at least the value of the observed r is "consistent with zero" , means that the most likely values 
of the inflationary density perturbation amplitudes, together with the CMB anisotropies induced 
by them, are infinitely large, or at least the data are "consistent" with such an infinity (see end 
of Sec lVIip . It goes without saying that we are not using this relationship. But we do use the 
relationship nt = — 1 which is a consequence of the superadiabatic generating mechanism, see 
Eq.dMl). 

The WMAP community seems to be satisfied with the data analysis which states that the CMB 
data can be described by a small number of parameters which include density perturbations, but 
with no necessity for gravitational waves. The usual logic in these derivations is first to find the 
best-fit parameters assuming that only the density perturbations are present, and then to claim 
that there is no much room left for inclusion of gravitational waves. This looks like being satisfied 
with the statement that most of what is known about the human race can be described by a small 
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number of parameters which includes one leg of individuals, but with no need for another. And 
when you propose to the data analysts that it is better to treat the data under the assumption 
that humans have two legs, they reply that this would be one extra parameter, and the proposer 
should be penalized for that. Anyway, in our analysis, we include the (inevitable) gravitational 
waves from the very beginning. 

The total TE signal is the sum of g.w. and d.p. contributions. Even if this sum is positive in the 
interval £ < 70, the effect of gravitational waves can (and expected to) be considerable. In this case, 
the amount of gravitational waves can be estimated through the analysis of all correlation functions 
together. However, if the total TE signal is negative in this interval of i, there is little doubt that a 
significant g.w. component is present and is responsible for the negative signal, because the mean 
value of the TE signal cannot be negative without gravitational waves (the issues of statistics are 
discussed below). 



It is intriguing that the WMAP team [3^ explicitly emphasizes the detection of a negative 
correlation (i.e. anticorrelation) at the multipoles near i ~ 30: "The detection of the TE anticor- 
relation near ^ w 30 is a fundamental measurement of the physics of the formation of cosmological 
perturbations...". The motivation for this statement is the continuing concern of CMB observers 
about the so-called 'defect' models of structure formation. These 'causal' models cannot produce 
any correlations, positive or negative, at ^ < 100. So, a detected correlation near £ 30 is an 
evidence againt them. At the same time, the better detected TE anticorrelation in the region of 
higher i « 150 could still be accomodated by the causal models, and therefore this more visible fea- 
ture is not a direct argument against these models. However, the motivation for the WMAP team's 
statement is not essential for the present discussion. Even if the available data are not sufficient to 
conclude with confidence that the excessive TE anticorrelation at lower £'s has been actually de- 
tected, it seems that the WMAP's published data (together with the published statement, quoted 
above) can serve at least as an indication that this is likely to be true. In the framework of the 
theory that we are discussing here, such an TE anticorrelation is a natural and expected feature 
due to relic gravitational waves, but of course this needs to be thoroughly investigated. 

The ensemble-averaged correlation functions do not answer all the questions. It is necessary to 
know what one would get with individual realizations of caused by individual realizations of 
random coefficients Cn, c^. The general unbiased estimator of Cf-^ is given by 

p,xx' _ ^ sr^ (x X'* ,„x*x'\ (oA\ 
^ ' m=-e 

where af^ depend in a complicated, but calculable, manner on mode functions, astrophysical input, 

49 



and, in general, on all random coefficients c^. (It was explicitely shown [48(] that the estimator 
DJ^ is not only unbiased, but also the best, i.e. the minimum- variance, estimator among quadratic 
estimators.) One can write, symbolically, 



-f-CXJ 

/ d'n J2 [fLin, s) c„ +/i^*(n, s) C] • (85) 



s=l,2 



Some features of the averaged XX' functions remain the same in any realization, i.e. for any 
choice of random coefficients c^, c^. One example is the absence of BB correlations for density 
perturbations. In the case of density perturbations, the BB correlation function vanishes not just 
on average, but in every realization. This happens because all the functions fgfn{n,s) in Eg. ([55]) 
are zeros. Another example is the retention, in every realization, of positive sign of the auto- 
correlation functions XX' , with X' = X. Indeed, expression (|85p may be arbitrarily complicated, 
but the estimator (j84|) for X' = X is always positive, as it is the sum of strictly positive terms. In 
different realizations the estimates will be different, but they will always be positive. That is, in 
any realization, including the actually observed one, the quantity Df-^ has the same sign as the 
quantity Cf-^. 

The important question is how often the TE estimations (j84p retain the same sign as the already 
calculated ensemble-averaged quantity Cj^. In other words, suppose a negative TE is actually 
observed in the interval i < 50. Can it be a statistical fluke of density perturbations alone, rather 
than a signature of gravitational waves ? In general, the answer is 'yes', but we will present 
arguments why in the problem under discussion the answer may be 'no', or 'very likely no'. 

To quantify the situation, it is instructive to start from two zero-mean Gaussian variables, say, 

and a^: 

{JJ) = a'^, {a^a^)=a'^, {Ja^) = pa'^, Q<p<l. (86) 

The probability density function for the product variable a^a^ involves the modified Bessel function 
Kq and is known as an exact expression 50']. In the situation, like ours, where p is close to 1, 
the probability of finding negative products a^a^ is small. A crude evaluation shows that, with 
probabilty 68% and higher, the values of a^a^ lie in the interval (0.3 — 1.7) around the positive 
mean value (a^a^), whereas for p = 1 the probability of finding a negative {a ^a^) is strictly zero. 
For 2i + l degrees of freedom the scatter around the mean value is expected to be much narrower. 
Nevertheless, if p 7^ 1, infrequent realizations with negative values of the product a^a^ are still 
possible. 
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We are probably quite fortunate in our concrete case of TE cross-correlation functions. It is 
true that af^ and af^ contain a large number of independent random coefficients Cn, multiplied 
by different deterministic functions fj^in^s), f^{n,s). However, in practice, these deterministic 
functions are such that the integrands in Eq. (j85p are quite sharply peaked at n ~ ^. Ideally, for 
every i, the random variables aj^, af^ become proportional to one and the same linear combination 
of random coefficients with |n| = n ^ i. The square of this combination is always positive, so the 
sign of the estimator DJ^ would be the same as the sign of the mean value Cj^ . This sign 
is determined by the known deterministic functions, not by statistics. Although this practical 
retention of the sign is not a strictly proven theorem, I think it is a very plausible conjecture. 

If the above-mentioned conjecture is correct, the negative sign of any observed TE correlation 
in the interval ^ < 50 is unlikely to be a statistical fluke of density perturbations alone, it should be 
a signature of presence of gravitational waves. One especially interesting interval of TE searches 
is the interval between £ ~ 30 and £ ~ 50. In this interval of i, the total signal Cj^ should be 
negative, whereas the contribution from density perturbations is still positive (see FiglTO] for a 
realistic example). The total TE'-signal is a factor 50, or so, larger than the expected i?i?-signal. 
All the logic described above suggests that if a negative TE is detected in this interval of ts it 
must have occured due to relic gravitational waves. 



XIII. PROSPECTS OF THE CURRENT AND FORTHCOMING OBSERVATIONS 



It is difficult to predict the future, but it seems to me that the Planck mission 
as the ground-based experiments, such as BICEP 
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53|, Clover 



5J, QUIET 



51[, as well 



551 ], have a very 



good chance of detecting relic gravitational waves through the TE anticorrelation described above. 
A level of the expected BB correlation is shown in Figj9j This level is natural for the discussed 
theory, but may be a little bit optimistic, as it is shown for a relatively high primordial spectral 
index n = 1.2. Observations with the help of space and ground facilities, and the measurement of 
all relevant correlation functions, should bring positive results given the expected sensitivities of 
those observations. It seems to me that the detection of relic gravitational waves in CMB is the 
matter of a few coming years, rather than a few decades. The analysis of the latest WMAP data 
provides serious indications of the presence of relic gravitational waves in the CMB anisotropies, 
see 0, 0- 



51 



Acknowledgments 



I am greatful to D. Baskaran and A. G. Polnarev for collaboration on the joint paper 
that was extensively quoted in this presentation, and to S. Weinberg for helpful comments on this 
manuscript. 



[1] C. W. Misner, K. S. Thorne, and J. A. Wheeler. Gravitation, Freeman and Co., San Francisco, 1973 
[2] L. P. Grishchuk. Zh. Eksp. Tear. Fiz. 67, 825 (1974) [Sov. Phys. JETP 40, 409 (1975)]; Ann. NY Acad. 

Sci. 302, 439 (1977); Pisma Zh. Eksp. Tear. Fiz. 23, 326 (1976) [JETP Lett. 23, 293 (1976)]; Uspekhi 

Fiz. Nauk 121, 629 (1977) [Sov. Phys. Usp. 20, 319 1977] 
[3] A. D. Sakharov. Zh. Eksp. Teor. Fiz. 49, 345 (1965) [Sov. Phys. JETP 22, 241 (1966)] 
[4] E. Schrodingcr. Physica 6, 899 (1939) 
[5] L. P. Grishchuk. Phys. Rev. 1)48, 5581 (1993) 

[6] L. Parker. Phys. Rev. Lett. 21, 562 (1968); Phys. Rev. 183, 1057 (1969) 

[7] B. Ya. Zeldovich. Impedance and parametric excitation of coupled oscillators, to be published in Physics- 
Uspekhi (2008) 

[8] L. D. Landau and E. M. Lifshitz. Mechanics. Oxford, Pergamon Press, 1975 

[9] P. Paradoksov. Uspekhi Fiz. Nauk 89, 707 (1966) [Sov. Phys. Usp. 9, 618 (1967)] 

[10] L. D. Landau and E. M. Lifshitz. The Classical Theory of Fields. Oxford, Pergamon Press, 1975 

[11] L. P. Grishchuk. Class. Quant. Gravity 10, 2449 (1993) 

[12] L. P. Grishchuk. Lecture Notes in Physics, Vol. 562: Gyros, Clocks, Interferometers: Testing Relativistic 

Gravity in Space, 167 (2001) 
[13] L. P. Grishchuk and M. Solokhin. Phys. Rev. D43, 2566 (1991) 
[14] S. Weinberg. The Quantum Theory of Fields. Cambridge U. Press, 1995 

[15] L. P. Grishchuk. Uspekhi Fiz. Nauk 175(12) 1289-1303 (2005) [Physics-Uspekhi, 48(12) 1235-1247 

(2005] , |arXiv:gr-qc/050401 8l 
[16] L. P. Grishchuk and A. G. Polnarev. In: General Relativity and Gravitation, 100 Years After the Birth 

of A. Einstein, ed. A. Held (Plenum Press, NY, 1980) Vol. 2, p. 393 
[17] D. Baskaran, L. P. Grishchuk, and A. G. Polnarev. Phys. Rev. D74, 083008 (2006) 
[18] Ya. B. Zeldovich and I. D. Novikov. The Structure and Evolution of the Universe. Chicago U. Press, 

1983 

[19] L. P. Grishchuk. Phys. Rev. D50, 7154 (1994) 

[20] L. P. Eisenhart. Continuous Groups of Transformations. Dover, 1933 
[21] J. M. Bardeen. Phys. Rev. D22, 1882 (1980) 

[22] J. M. Bardeen, P. J. Steinhardt, and M. S. Turner. Phys. Rev. D28, 679 (1983) 



52 



[23] V. N. Lukash. Pis'ma Zh. Eksp. Tear. Fiz. 31, 631 (1980) [JETP Lett. 6, 596 (1980)] 
[24] G. Chibisov and V. Mukhanov. Mon. Not. R. Astr. Soc. 200, 535 (1982) 
[25] M. Sasaki. Prog. Theor. Phys. 76, 1036 (1986) 



[26] W. Unruh. Cosmological long wavelength perturbations., arXiv:astro-ph/9802323 
[27] S. Weinberg. Phys. Rev. D72, 043514 (2005) 



[28] R. Flauger and S. Weinberg. Phys. Rev. D75, 123505 (2007) ( |arXiv:astro-ph/0703179| ) 
[29] S. Weinberg. Phys. Rev. D69, 023503 (2004) 

[30] V. N. Lukash. On the relation between tensor and scalar perturbation modes in Friedmann cosmology, 
|arXiv:astro-ph/0610312| 

[31] D. N. Spergel et al. WMAP Collaboration, Astrophys. J. Suppl. 170, 377 (2007) 
( |arXiv:astro-ph70603449| ) 

[32] L. Page et al. WMAP Collaboration, Astrophys. J. Suppl. 170, 335 (2007) f arXiv:astro-ph0603450"| ) 
[33] L. Pagano, A. Cooray, A. Melchiori, and M. Kamionkowski. Red Density Perturbations and Inflationary 

Gravitational Waves, larXiv:0707. 25601 (In this fresh paper one can also find some previous history of 

derivation and appUcation of the infiationary 'consistency relation'.) 
[34] H. V. Peiris et al. Astrophys. J. Suppl. 148, 213 (2003) 

[35] S. Weinberg, website: www.ainacad.org/publications/bulletin/summer2003/weinberg.pdf 
[36] S. Bose and L.P. Grishchuk. Phys. Rev. D66, 043529 (2002). 
[37] R. K. Sachs and A. M. Wolfe. Astrophys. J. 147, 73 (1967). 
[38] G. Dautcourt Month. Not. R. Astr. Soc. 144, 255 (1969) 

[39] L. P. Grishchuk and Ya. B. Zeldovich. Astron. Zh. 55, 209 (1978) [Sov. Astron. 22, 125 (1978)]. 
[40] A. G. Polnarev. Sov. Astron. 29, 607 (1985) 

[41] S. Chandrasekhar. Radiative Transfer. New York, Dover Publications Inc., 1960 

[42] M. Zaldarriaga and U. Seljak. Phys. Rev. D 55, 1830 (1997) 

[43] M. Kamionkowski, A. Kosowsky, and A. Stcbbins. Phys. Rev. D 55, 7368 (1997) 

[44] M. M. Basko and A. G. Polnarev. Sov. Astron. 24, 268 (1980) 

[45] U. Seljak. Astrophys. J. 482, 6 (1997) 

[46] W. Hu and M. White. New Astron. 2, 323 (1997) 

[47] A. M. Yaglom. An Introduction to the Theory of Stationary Random Functions. Prentice-Hall, Engle- 

wood Cliffs, NJ, 1962 
[48] L. P. Grishchuk and J. Martin. Phys. Rev. D56, 1924 (1997) 
[49] R. Crittenden, D. Coulson, and N. Turok. Phys. Rev. D52, R5402 (1995) 
[50] L. P. Grishchuk. Phys. Rev. D53, 6784 (1996) 



[51] The Planck Collaboration. arXiv:astro-ph/0604069 



[52] B. Keating and N. Miller. New Astron. Review, 50, 184 (2006) 

[53] J. M. Kovac and D. Barkats CMB from the South Pole: Past, Prsent, and Future, [arXiv:0707.1075l 
[54] A. Taylor et al. In: Proceedings of 39th Recontres de Moriond (Frontieres, 2004) 



53 



[55] C. Lawrence. The QUIET Experiment, (http://hea.iki.rssi.ru/Z-90/). Invited talk at the conference 

"Zeldovich-90" , Moscow, 2004 
[56] W. Zhao, D. Baskaran, and L. P. Grishchuk. Phys. Rev. D79, 023002 (2009) 
[57] W. Zhao, D. Baskaran, and L. P. Grishchuk. Phys. Rev. D80, 083005 (2009) 

[58] All inflationary quotations in the present article are taken from real publications. However, I do not 
see much sense in giving precise references. I am confronting here something like a culture, distributed 
over hundreds of publications, rather than a persistent confusion of a few authors. This situation was 
once qualified 26| as "the controversy between Grishchuk and the rest of the community" . 



54 



